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Results from general relativistic theoretical Monte Carlo computer simulations of Compton scat- 
tering and e~e + pair production processes in the ergosphere of a supermassive (~ 10 s Mq) rotating 
black hole are presented. Particles from an accretion disk fall into the ergosphere and scatter off 
particles that are in bound equatorially and nonequatorially confined orbits. The Penrose mecha- 
nism, in general, allows rotational energy of a Kerr black hole to be extracted by scattered particles 
escaping from the ergosphere to infinity (i.e., large distances from the black hole). The results of 
these model calculations show that a form of the Penrose mechanism is capable of producing the 
observed high energy particles (up to ~ GeV) emitted by quasars and other active galactic nuclei 
(AGNs), without the necessity of the external electromagnetic field of the accretion disk. Impor- 
tantly, these model calculations show that the Lense-Thirring effect, i.e., the dragging of inertial 
frames into rotation, caused by the angular momentum of the rotating black hole, results in a grav- 
itomagnetic (GM) force being exerted on the scattered escaping particles. Inside the ergosphere, 
where this dragging is severe, in appears that the GM vector field lines are frame dragged into the 
„Cj ■ positive azimuthal direction, i.e., the direction of rotation of the black hole. The resulting GM force 

acting on the Penrose scattered particles produces symmetrical and asymmetrical (or one-sided) 
particle emissions in the polar direction, consistent with the astrophysical jets observed in radio 
00 ' strong AGNs. Note, these Penrose processes can apply to any size rotating black hole. 
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(N '. I. INTRODUCTION 

In 1918, it was discovered in the equations of general relativity by Thirring and Lense [1] that the angular momentum 
psj I of a massive body causes inertial frames to be dragged around in the direction that the mass is rotating. In 1969, 
Roger Penrose [2] noted that, inside the ergosphere of a rotating black hole, this frame dragging is so extreme that it 
can provide a way to extract rotational energy (we now call this the Penrose mechanism). In this paper, I present the 
results of the first complete, fully general relativistic, four-dimensional, theoretical and numerical model calculation 
using the Penrose mechanism to extract high energy electrons and photons from a rotating black hole. Penrose 
scattering processes in the ergosphere allow high energy particles to escape to infinity (far away from the black hole) 
with a portion of the rotational energy of the Kerr black hole (KBH), and orbital energy- momentum produced by 
its strong gravitational field. This model calculation not only confirms the results of energy extracted from Penrose 
• • ■ Compton scattering processes given in [3,4], however, extracting higher energies when nonequatorially confined targets 
are introduced, but this model gives space momenta of individual scattered particles, thus, in this sense, making it 
a complete calculation. Importantly, the Penrose pair production, by photons in bound unstable trapped orbits with 
photons on radially infalling geodesies, at the photon orbit, allow electron-positron (e~e + ) pairs to escape with energies 
as high as ~ 4 GeV (implying a Lorentz factor 7 e ~ 10 4 ), or higher depending on properties of the accretion disk. Up 
until the investigation reported here and the related paper of [5], such Penrose pair production (PPP) at the photon 
orbit had yet to be proposed. The reason for this is probably because the trajectories for massless and material 
particles, in nonequatorially confined orbits, i.e., orbits not confined to the equatorial plane of the KBH, had to be 
solved first to obtain analytical expressions for the conserved energy E and conserved azimuthal angular momentum 
L, as measured by an observer at infinity; this is done in [5], and the results are presented again here in Section fl A| . 

The classical Penrose process utilizes the existence of retrograde particle orbits (with respect to the rotation of 
the KBH) in the ergosphere, for which the energy, as would be measured by an observer at infinity, is negative [2,3]. 
Such orbits do not come out to infinity. However, it is possible for a particle, say p±, that has fallen inwardly from 
infinity into the ergosphere to scatter off another particle, say P2, initially in a direct orbit inside the ergosphere. If 
the orbit of the scattering particle P2 changes into a retrograde orbit (of negative energy), then the scattered particle 
Pi can escape to infinity with more mass-energy than the sum of the initial energies of p\ and pi- Since the orbit of 
the initially bound particle P2 is dependent on the angular momentum of the KBH and the curvature of spacetime 
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defined by the mass of the hole, when pi gives up energy to the escaping particle p± and falls into the event horizon, 
the KBH loses energy in the form of rotational energy. 

One plausible class of Penrose processes occurs when particles already inside the ergosphere (say particles of an 
accretion disk) undergo local relativistic scatterings [6,7]. Such types are considered in this investigation. If one of the 
scattering particles is initially in a bound orbit, then it is possible for the other initially unbound scattered particle, or 
the new particles created in the scattering process, to escape to infinity with rotational energy-momentum from the 
KBH, extracted either directly (classical Penrose process) or indirectly (quasi-Penrose process), as defined in [5]. This 
process allows scattered ergospheric particles to (1) escape to infinity with more mass-energy than they would have 
had if the scattering occurred outside the ergosphere and (2) escape to infinity with mass-energy initially trapped by 
the KBH, mass-energy near the event horizon that had no other way of escaping, except by such processes as these 
Penrose processes; otherwise this trapped mass-energy would eventually be accreted into the black hole. That is, 
since, in general, nothing can come out of the black hole (including for orbits near the event horizon), then a physical 
process near the hole, such as a Penrose process, is necessary to inject a trapped or "plunging" orbiting particle into 
an escaping orbit [7]. In other words, these so-called Penrose processes provide a way for blueshifted mass-energy, 
initially trapped by the KBH, to escape to infinity. 

The model calculations discussed in this paper consist of using the Monte Carlo method to treat Penrose processes 
in the ergosphere of a KBH. The processes investigated are Penrose Compton scattering (PCS), 7-ray — proton PPP 
(7P — > e ~ e+ p)i an d 7-ray — 7-ray PPP (77 — > e~e + ). These three processes are assumed to occur in material that 
has fallen into the ergospheric region of a rotating black hole from a surrounding accretion disk. The Monte Carlo 
method is applied to the cross section of each scattering event to determine values of the scattering angles and final 
energies as measured by local observers. 

In [5], I present the detailed model calculations for the above Penrose scattering processes; the final outcome is 
the four-momenta of the scattered particles as measured by an observer at infinity. The important properties and 
relations found from the model calculations of [5] are described in this paper. Now, of the three Penrose scattering 
processes investigated, the most rewarding ones were the PCS and the PPP (77 — > e~e + ): these two processes are 
specifically discussed in this paper. Unfortunately, the PPP (7P — > e~e + p) as suggested in [8,9] did not allow any 
of the scattered e~e + pairs to escape, mainly because of the large inward radial momentum acquired by the assumed 
radially infalling incident photon at the scattering radius, too large for the pairs to be scattered outward by the target 
orbiting proton: I refer the reader to [5] for details concerning this scattering process. 

The purpose of this paper is to discuss some of the astrophysically interesting features found by looking carefully at 
the emission spectra of the Penrose scattered particles. Two of the most important features of the scattered particles 
are the following: (1) large fluxes of relativistic e~e + pairs with energies as high as ^4 GeV can escape; and (2) the 
rotation of the KBH can cause the scattered particles to escape with asymmetrical distributions, above and below 
the equatorial plane, thus, naturally aiding in the production of one-sided jets of relativistic particles. Feature (2) 
is a general relativistic effect caused by the gravitomagnetic field [10] (i.e., the gravitational analog or resemblance 
of a magnetic field), due to the angular momentum of the rotating black hole. This gravitomagnetic (GM) field, as 
described in Section [ID , causes the dragging of inertial frames, while at the same time produces a force which acts on 
the momentum vector of a particle in its field, and which, for small velocities c), is linear in the particle's velocity. 
Both of the above features are important ingredients in what would be needed to explain the observed spectra of 
quasars and other AGNs. Note, feature (2) implies that, near the event horizon the expected reflection symmetry of 
the Kerr metric, above and below the equatorial plane, appears to be broken due to the GM force field. 

It has been suggested [11] that the one-sided or asymmetrical jets of radio strong quasars and radio galaxies may 
be intrinsically related to the energy source. The reason for this suggestion is that many of the observed one-sided jet 
distributions of AGNs do not conform with that predicted if the one-sidedness is due solely to relativistic beaming — of 
the emission axis of the jet into a direction near the line of sight of the observer — wherein Doppler boosting of the 
radiation causes the approaching jet to be apparently brighter than the receding jet. Yet, up until the investigation 
reported here, no model has been proposed for the energy source that would exhibit such asymmetrical properties. 
Herein, however, from these Penrose scattering processes [5], I present evidence that the rotating black hole core 
energy source may be producing two of the important observed features of AGNs: the large energy fluxes of e~e + 
pairs inferred to emerge from their cores, and asymmetry in their jet and counter jet distributions. Moreover, findings 
by Dennett-Thorpe et al. [12], after studying the asymmetry of jets, lobe length, and spectral index in a sample of 
quasars, with well defined jets, led them to conclude that the simplest models of relativistic beaming fail to explain 
the "present observations"; and that only if the intrinsic spectrum is curved (i.e., with the spectral index increasing 
with frequency), such that the approaching lobe is seen at a significant lower frequency, will asymmetry of the right 
sense be produced. This somewhat supports the findings presented here, that asymmetry in the jet and counter jet 
is intrinsic to the energy source, suggesting that, relativistic beaming may just serve as an important enhancement 
mechanism. Further, recent observations of so-called microquasars in our Galaxy show some asymmetry in the jet 
and counter jet distributions [13]; and not all of the observed asymmetry can be explained by relativistic beaming, 
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suggesting that some asymmetry is intrinsic [14]. 

The remaining structure of this paper is as follows: In Section 0, 1 present general formalisms and descriptions of the 
Kerr metric, particle orbits, the Penrose Compton scattering (PCS), the Penrose 7-ray — 7-ray pair production [PPP 



(77 — * e ~ e +)]i an d the gravitomagnctic field (GM). Presented in Section III is a discussion of the gravitomagnetic 
field as it relates specifically to PCS and PPP (77 — ► e~e+). Finally, in Section I conclude with an overview, 
including general discussions of how well these model results agree with astronomical observations, and the role of an 
external accretion disk magnetic field. Also, suggestions are made for further investigations. 

II. FORMALISM 

A. The Kerr Metric and Particle Orbits 

For completeness, review, and to be referred to in the following sections, the spacetime separation between events 
near the KBH is defined by the Kerr metric [15]. The Kerr metric written in Boyer-Lindquist coordinates [16] and in 
geometrical units (G = c = 1) is 

ds 2 = -e 2v dt 2 +e 2 ^{d$- 

with 

e 2 » = 

e 2V = 

e 2 Ml = 
e 2„2 = 
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a is the angular momentum per unit mass parameter, M is the mass of the black hole, and tu gives the frame dragging 
velocity. The quantities A and E are defined by 

A = r 2 - 2Mr + a 2 (8) 

and 

E = r 2 + a 2 cos 2 6, (9) 

respectively. In general, the parameter a can have values < (a/M) < 1, values which allow for the existence of 
an event horizon. For a KBH, the event horizon is located at r = r + = M + I M 2 — a 2 ) , the larger root of the 



equation A = 0. Here, in these Penrose processes, a canonical KBH is used, with its limiting value a/M = 0.998, as 
defined in [17], in the investigation of an accretion disk around a KBH. It is found that, if a/M is initially very close 
to 1 (0.999 < a/M < 1), a small amount of accretion (AM/M < 0.05) through a disk quickly spins the hole down 
to a limiting state a/M ~ 0.998. Conversely, if a/M is initially below this limiting value, accretion spins the hole up 
toward it. 

Upon approaching a KBH from infinity, a limit is reached where the angular momentum of the KBH causes inertial 
frames to be dragged around in the direction that the black hole rotates. This limit is characterized by the vanishing 
of the gtt component in Eq. ([!]): 

ff« = -fl--g-J=0. (10) 
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More precisely, this limit is the larger root of g t t, given analytically by 
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The region between this limit and the event horizon is called the ergosphere. 

Inside the ergosphere, the Kerr metric in the Boyer-Lindquist coordinate frame (BLF), i.e., frame of the observer at 
infinity, does not allow an observer to be stationary (in the sense of the observer being at rest with r, O, $ = constant) 
because of the dragging of inertial frames. For this reason, physical processes are difficult to describe in the BLF. 
In order to examine physical processes inside the ergosphere, Bardeen, Press, and Teukolsky [7] devised a frame of 
reference called the local nonrotating frame (LNRF). Observers in this frame rotate with the KBH in such a way that 
the frame dragging effect of the rotating black hole is canceled as much as possible. In the LNRF, special relativity 
applies since locally spacetime has Lorentz geometry. The transformation laws, for the covariant components of a 
four-momentum, between the BLF and the LNRF tetrads are given in [3,5]. Note, in the LNRF, g^ v = r}^ = 
Minkowski metric components. 

However, it is better to use the BLF when describing the general orbits of particles (including photons) near the 
event horizon of the KBH. The BLF admits three constants of motion as measured by an observer at infinity [18]. 
In terms of the covariant components of the particle's four-momentum [P M = (P r , P e , P$, P t ) ; = r, 0,$,t)l at 
some instant of time, the conserved quantities are E = —P t = total energy, L = P$ = component of the angular 
momentum parallel to the symmetry axis, and 



Q = Pi + cos 2 6 



Mo 



E' 



L 2 



(12) 



where fj, is the rest mass energy of the particle, which is a trivial fourth constant of motion. (Q is sometimes referred 
to as Carter's constant.) The value of Q is zero for particles whose motions are confined to the equatorial plane. 
The nonzero values of Q belong to particles which are moving in the ±e© direction and/or are not confined to the 
equatorial plane. 

The general expressions for the conserved energy E and azimuthal angular momentum L of massless (fi Q — 0) and 
material particles on orbits not confined to the equatorial plane (i.e., with Q > 0) in the curved spacetime of a KBH 
are given below. For direct orbits of constant radius r [5] 
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Note that, when Q = 0, Eqs. ( |l3| ) and ( |l4| ) reduce to the equations for equatorially confined orbits given in [7,5]. 
The quantity y/Q (= Pq) versus the conserved energy E [Eq. (|l3|)] is plotted in Fig. 1(a) for the photon orbit at the 
radius r p h, and in Fig. 1(b) for the electron orbits at the ma rgina lly boun d and the marginally stable radii, r m b and 
r ms , respectively (these radii are discussed further in Sections II B and II C ). Importantly, the derivations of Eqs. (|T^ ) 
and ([Ml, derived in [5], now allow for Penrose processes be investigated in a practical way. That is, the motion of 
particles, above or below the equatorial plane (in the ±e© direction), can now be considered for the target particles 
as well as for the scattered particles, with E and L given separately by Eqs. (|l^) and (|lj), respectively, wherein Q is 
the independent variable for a specific black hole of mass M and angular momentum parameter a, at constant radius 
r. 

When Q > 0, Eqs. ( |l3| ) and Jl4| ) describe nonequatorially confined "spherical-like" orbits [19]. One of these 
spherical-like orbits, near the event horizon, consists of a particle repeatedly passing through the equatorial plane 
while tracing out a helical belt lying on a sphere at constant radius. The belt width between the maximum and 
minimum latitudes that the particle achieves, in general, increases with increasing Q, where Q is given by Eq. (|l^). In 
these spherical-like orbits, which in the limit of large radius goes asymptotically to a Keplerian circle, as the orbital 
radius approaches that of the event horizon, the dragging of the line of nodes, in the sense of the spin of the black hole, 
increases without limit. (The nodes are points at which the orbit, in going between negative and positive latitudes, 
intersects the equatorial plane.) Inside the ergosphere, during the time that a particle makes one oscillation in latitude, 
it will be swept through many complete azimuthal revolutions. Consequently, an orbit near the event horizon will 
have a helix-like shape with the axis parallel to the spin of the KBH. In the Penrose processes investigated here, it is 
assumed that a collision takes place, for a given observer, when the target particle with Q > passes through one of 
the nodes. 

After the scattering events, not all of the particles escape to infinity. There are certain conditions that a particle 
must satisfy in order for it to be on an outgoing orbit. A thorough discussion of the escape conditions for material and 
massless particles is given in [5] . Generally speaking, these escape conditions compare the four-momentum components 
E', L', and Pq (= y/Q 7 ) of the scattered particle to the values E 7 L, and Q of a possible turning point, to see if the 
particle is on an escaping orbit. The expressions for E and L are given by Eqs. ( |l3| ) and (0). Note that, for the 
nonequatorially confined (Q > 0) spherical-like orbits, bound particles can have energies E//i > 1, because of an 
additional escape condition that must be satisfied. 



B. Penrose Compton Scattering (PCS) 



In the PCS process, a photon of initial energy E p h, as measured by an observer at infinity (i.e., in the BLF), is 
assumed to be emitted inside the ergosphere — emitted from the innermost region of the accretion disk, and follows a 
null geodesic trajectory before colliding with an equatorial or a nonequatorially confined orbiting electron, located at 
a specific radius of energy E e and azimuthal momentum L e [= (P e )*] given by Eqs. ( |l3| ) and (pT[), respectively. The 
motions of the incoming or incident photons considered in these scattering processes are those moving radially inward 
along the equatorial plane, as in [3]. Now, some of the photons, after being scattered by the electrons, eventually 
escape to infinity. For the scattered photons that are allowed to escape to infinity, an observer there would observe a 
low-energy (in most cases < /z e , where \x e is the electron rest mass energy) photon being scattered by a direct orbiting 
electron, after which the photon comes out with a higher energy (inverse Compton scattering). Subsequently, the 
target electron may recoil to another direct orbit of lesser energy, or the electron may be put on a retrograde orbit of 
negative energy. In both cases, the target electron gives up energy as measured by an observer at infinity. However, to 
a particular observer in the LNRF this is just a normal Compton scattering process in which the photon loses energy 
to the electron, since the photon arrives at this local frame with initial energy higher than ~ /i e [5] . The infall of the 
final negative energy electron results in an observer at infinity measuring a decrease in the rotational energy of the 
KBH. This inverse Compton scattering process is different from its flat spacetime counterpart. In the flat spacetime 
process, cold photons are heated by hot electrons to the temperature of the electrons. In the Penrose process, which 
occurs in the ergosphere, generally, photons are heated by the rotational energy of the black hole and the curvature 
of spacetime [3]. 

The conditions that the infalling photon encounters in its free fall through the ergosphere before arriving at the 
designated scattering radius are ignored in these calculations. Otherwise, conditions such as the particle density in the 
ergosphere and radiative transfer effects would have to be incorporated into the calculations along the null geodesic 
of the photon. It is assumed that these conditions cause little qualitative change in the results. The validity of this 
assumption can be found in [4], where it is shown that the variation in the density distribution does little to change 
the energies of the escaping photons. The photon arrives at the radius in which the scattering take place with an 
initial covariant four-momentum (P p h)p = [(Pph)n (Pph)e> -^ph, — Pph], as measured by an observer at infinity. For a 
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photon falling radially inward along the equatorial plane, (P p h)e — -^ph = 0. 

The initial energies -Eph, used in the model calculations presented here, are monochromatic for a given distribution 
of 2000 incident infalling photons. One scattering event per photon is assumed. The energies of the incident photons 
(0.511 MeV— 1 MeV) are consistent with the bistable two-temperature accretion disk [8,20,21], referred to as the thin 
disk/ion corona accretion model (see [5] for further details concerning the disk model and the use of thermal initial 
particles). 

The spatial distribution of the target electrons is assumed to be specified as circular orbiting rings of completely 
ionized plasma. Two radii are considered for the model electron rings; r ms (— 1.2A^f), the radius of the marginally 
stable circular orbit and r m b (— 1.089M), the radius of the marginally bound circular orbit, for a canonical KBH 
(a/M = 0.998). Note that, the outer ergosphere boundary, the photon orbit, and the event horizon of a canonical 
KBH are located at r = 2M, r p h ~ 1.074M, and r + ~ 1.063M, respectively, in the equatorial plane. The region 
between the radii r ms and r m b represents the innermost possible periastrons (radial turning points) for bound unstable 
orbiting material particles. If given a sufficiently large inward perturbation, a particle in this region will eventually 
fall through the event horizon. Moreover, r m b is the minimum periastron of all bound orbits of material particles, 
before plunging directly into the black hole. 

After some manipulations of the results given in [5], eq. (3.47), the final four-momentum components of the PCS 
photon as measured by a BLF observer, can be expressed by the following: 
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where e u , 



, e Ml , e M2 , and w are given by Eqs. (|2|) through (J6J); the other parameters will be defined in the following 
paragraph. These four-momentum components result from scattering off equatorial target electrons (i.e., targets 
confined to the equatorial plane). See [5] to derive the general ex pre ssions for both equatorially and nonequatorially 
confined targets. Note that, the arccosine term in Eqs. (|l6| ) and ( jl7fe defines an angle between and 2tt. 

The other parameters in the above four- momentum components [Eqs. ( |l5| ) through (18)] are defined as follows [5]: 
For equatorial orbiting target electrons, j3 e = (f3 r , (3q, /?$) = [0,0, (w e )*], and 7 e = (1 — 0l)~ x ^ 2 , with 



(v e )<s> = (3$ = 



E„ — ujLf. ' 



(19) 



where E e and L e are the conserved energy and azimuthal angular momentum as measured by an observer at infinity, 
given by Eqs. (|l^) and (|l4|), respectively, of course with Q e — for equatorially confined targets (Q e ^ for 
nonequatorially confined targets); (u e )$ is the orbital velocity of the electron relative to the LNRF. Note that, in [5], 
I refer to the lab frame (LF) as the frame of a general observer in the LNRF; here I simply refer to it as the LNRF. 
The final energy of the PCS photon £p h , as measured by an observer in the electron rest frame (ERF), is given by 
the familiar expression for Compton scattering in flat spacetime (see [5]): 

R 



cos S R ) 
j e e~"E ph 



1 + (7 e e-" E ph / /z e )(l- cos S R )' 
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where e^ h is the initial energy of the photon in the ERF; S R is the angle between the initial and final space momentum 
vectors of the photon in the ERF. The ERF and the LNRF are related by a Lorentz transformation, with relative 
frame velocity (u e )$- The initial and final polar angles, 9 R h and 9 R h , of the incoming and scattered photon space 
vectors, respectively, in the ERF [in which the coordinate system is centered on the electron at rest, with the polar 
axis (+e z ) of this Lorentz frame pointed in the positive e$ direction of the BLF; see eqs. (3.3), (3.5), (3.6), (3.21), 
and (3.24) of [5]], are found to be 



7 ph 



arctan 



Oft)? + (*£)! 



(Pph)* 



1/2 



tt — arctan 



= 7r — arctan 



2 mlMar sin 2 6 



ujLl(r 2 + a 2 cos 2 8) 



1/2 



(21) 



using (p* ) r = -e-'Bph, (p R h ) e = 0, (p« )* = -^fatT"^ and eqs. (2.3), (2.10b), (2.12), (3.15) from [5]]; and 



9 R h = arccos(cos 9 R h cos (5 fi + sin 5 R sin cos a ) 



(22) 



[see also the azimuthal angle 0^ [eq. (3.39)] of [5], given again here as Eq. Ji4|)], where (Pp h ) rj (Pph)©> anc ^ (Pph)* are 
the initial space momentum components of the photon in the ERF; note, the tt term has been properly added since 
(Pph)* < 0. The angles 9 R h , (p R h define the direction of the space momentum vector of the scattered photon in the 
ERF (see eq. (3.24) of [5]). The angles S R and a R are the polar and azimuthal angles of the PCS photon relative to 
the direction of the initial photon: these angles are found by applying the Monte Carlo method to the Klein-Nishina 
cross section in the E RF (as done in [5] ) . Equations (^l[) and ( p2|) will be important in the discussions of the CM 
field, in Sections II D and III, because of the frame dragging term u contained in 9 R h [cf. Eqs. ( Jl5| ) through (|l8|), (|lS|), 
©], and ©]. 

In [5], I present numerous cases of the obtained four- momentum spectra of the PCS photons. Below in Section III, 
however, I present only spectra that are relevant to the discussion at hand — i.e., of how the GM field affects the distri- 
butions of the scattered photons. Note that, the conserved energies E e and angular momenta L e of the nonequatorially 
confined target electrons, assumed in these calculations, are consistent with the electrons produced in Eilek's hot ac- 
cretion disk model [8,20], described in details in [5]. The Q e values used in Eqs. ( |l3| ) and (14j) for the nonequatorially 
confined target electrons correspond to orbits with polar latitudinal angles ^30°, above and below the equatorial 
plane, consistent with what would be expected in thin disk/ion corona (or torus) accretion models. 



C. Penrose Gamma Ray-Gamma Ray Pair Production (PPP) 

This PPP (77 — > e~e + ) process consists of collisions inside the ergosphere between radially infalling photons [5] 
and circularly orbiting photons that are bound at the radius of the photon orbit [7]: 

r ph = 2M{1 + cos [(2/3) arccos (-o/M)]}, (23) 

where r p h represents an unstable circular orbit. On such spherical-like orbits, a photon circles indefinitely; however, 
a small perturbation can cause the photon to either escape to infinity or fall into the black hole. This orbit is the 
innermost boundary of circular orbits for all particles. The initial energies assumed in these calculations as measured 
by an observer at infinity, for the infalling photon and the orbiting photon, Eyi and E 7 2> respectively, should be in the 
following ranges to be consistent with energies produced by thin disk/ion corona models, and with energies satisfying 
the threshold energy conditions for these reactions to occur [5]: 

3.5 keV < E 7 i < 50 kcV and 3.4 MeV < E l2 < 4 GeV, 

where E 7 2 is given by Eq. ( |l3| ) for the appropriate chosen value of Q 7 2. Note that, the target photons, E^, have 
orbits not confined to the equatorial plane. These nonequatorially confined orbits allow for the photon's energy at 
the photon orbit to be used in a feasible way, whereas, for equatorially confined orbits, Eqs. (|l3|) and (|l4|), for E and 
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L, respectively, equal zero [cf. Fig. 1(a)]. Or expressed in conventional terms, E and L, when Q=0, define the orbit 
of infinite energy and infinite azimuthal angular momentum per unit rest mass energy, respectively [7,5]. 

The high energy range for E l2 is chosen based on physical processes found to occur in the two-temperature bistable 
thin disk/ion corona accretion flow. More recently such type ion corona accretion flows have been called advection 
dominated accretion flows (ADAFs). ADAFs were predicted in Lightman's speculation of what might take place 
as a result of thermal and density instabilities [22]. These instabilities can cause the thin disk to form high and 
low density rings in the inner region and an increase in temperature from ~ 10 5-7 K to ~ 10 9-12 K [8,20,21,22]. 
Astrophysicists over the years, since Lightman's time-dependent disk model calculation of the early 1970s [22], and 
subsequent models [8,20], have continued to find properties of ADAFs [23]. The most important is the proton-proton 
collisions in which neutral pions ir° are emitted [8,20,23]. The neutral pion ir° subsequently decays into two energetic 
7-rays with energies up to ~ 100 MeV if the protons are assumed to have a blackbody distribution, or higher for an 
assumed power-law distribution [23]. Such photons could populate the photon orbit as follows. These inf ailing photons 
would be blueshifted in energy-momentum at the photon orbit by a factor ~ e~ v ~ 52 (due to the rotating black 
hole producing the Lense-Thirring effect), where e~ v is the so-called blueshift factor given by Eq. (g) [cf. eq. (2.8d) 
of [5]]. This means that high energies up to ~ 5 GeV (or higher) can possibly exist in the ergosphere at or near 
the photon orbit. However, it requires a process, such as the Penrose mechanism, to inject these highly blueshifted 
trapped orbiting photons into escaping orbits [7] — in this PPP (77 — ► e~e + ), they escape in the form of energetic 
e~e + pairs. Now, it is well known that for rotating black holes with a/M — 0.998, the accretion disk extends inward 
to ~ r ms (~ 1.2M), inside of the radius of the ergosphere (r Q = 2M). So, the 7 rays produced in Eilek's hot accretion 
disk model [8,20], from the decay of neutral pions ir° — ► 77, described in details in [5], with energies peaked around 
~ 75 MeV, could very well be seed photons for these 77 — > e~e + reactions, photons that have become blueshifted 
and bound at the photon orbit [24]. In addition, the PCS processes of [5] show that some of the scattered photons can 
also be seed particles for these reactions. Further, these photon collisions can occur at other radii besides the photon 
orbit, at, say another turning point, thereby increasing the probability of a collision. Again, it is important to point 
out that once the photons have been blueshifted to energies given by £^2, and become bound at the photon orbit, 
they have no other way of escaping, except by some perturbing physical interaction, such as these PPP (77 — > e~e + ) 
processes. Note that, in [5], these PPP (77 — ► e~e + ) processes are defined as quasi-Penrose processes since they are 
somewhat different from classical Penrose processes — which utilize the retrograde orbits of particles in the ergosphere. 

Before proceeding, I want to discuss the assumptions made above and their validities. First, it is assumed that the 
blueshift in the energy of an infalling orbiting photon as measured by a LNRF observer is approximately equal to that 
measured by a BLF observer. This is a valid assumption because, for a photon in a circular bound orbit about the 
KBH, there is little difference between the energy measured locally in the LNRF [eq. (2.8d) of [5]] and that measured 
by an observer at infinity (i.e., in the BLF) [Eq. (jl3[)]; for example, at the photon orbit r p h, the ratio of the energy 
measured by an observer in the LNRF, e* h , to that measured by an observer in BLF, E* h , is ~ 1.055. Second, it 
is assumed that the photons populating the photon orbit have been created by some prior processes (disk instability 
and/or Penrose), which yielded them with appropriate energies and momenta. We know from [5] that the condition 

, » Qph Q^h 

^ph h ph 

guarantees the existence of a turning point outside the event horizon for which the radial motion goes to zero and 
changes direction for photons emitted not confined to the equatorial plane [i.e., (P ph )e 7^ 0]. The asterisks indicate 
the nonequatorially confined circular photon orbit, where r = r p h [see Eq. (p3|)], which is a potential turning point, 
and the primes indicate a photon resulting from some prior process, say PCS. Once the condition above has been 
met, indicating the existence of a turning point outside the horizon, to find out whether or not the photon orbit is 
a true turning point, for a particular photon, it seems plausible to set Q ph = Q ph , in the the condition above, and 
let E' ph — ► e~ v E' ^ = E p ^, as the photon spirals inward toward the photon orbit, obtaining the requirement that 
E ph < E* h for the existence of a turning point at r = r p h [where e~ v is the gravitational blueshift factor, again, 
given by Eq. (^) , due to the local BLF orbital energy gain; cf. also Eq. ([[])] . After the photon has reached its turning 
point, it seems reasonable that it escapes outward with "redshifted" energy given by E' ph = e u E ph as measured by an 
observer at infinity, i.e., in the BLF. Although at this point we do not know the exact detailed path of the photon to 
and from its turning point, but, we can theorize that it looses its acquired inspiraling local orbital energy on its way 
back outward, before escaping to infinity, as indicated in the above relationship between E' ph and -EpV Further, notice 
that the condition above is independent of L' ph . The other condition that independently guarantees the existence 
of a turning point in P ri however, for photons confined to escape along the equatorial plane, as well as for those 
escaping not confined to the equatorial plane, involves L' ph ; see [5] for details and a complete description of the escape 
conditions. In [5], it is found that many of the PCS photons have appropriate E' h and Q' h values to populate the 
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turning point trajectories at or near the photon orbit [cf. Figs. 3 and 4 to Fig. 1(a)]. Because of this, it is not 
unreasonable to expect the photons produced by the ir° — > 77 decays (discussed above) to acquire appropriate E' ph 
and Q' ph values, particularly, if the process, starting from the proton-proton scattering, is properly considered as a 
Penrose process inside the ergosphere. Note, the radii of other turning points can be found from Eqs. @ and @), 
by satisfying the condition E ph = e~ u E' ph < E ph , for the specified "iso-Q" value orbit ( = circular orbit of equal Q), 
where now the radius of the turning point can exist anywhere near the event horizon, say between r m b > r > r p h for 
the inward PCS photons. This means that such photons satisfying a condition for the existence of a turning point can 
have energies at these turning points corresponding to blueshift factors in the range of 32 < e~ v < 52. The largest 
factors, however, will be for photons nearest the photon orbit. In these calculations, the investigation centers around 
those photons with turning points at the photon orbit where the highest energy can be extracted. Yet, consideration of 
all of the possible turning points will increase the probability of a PPP (77 — ► e~e + ) process occurring, as mentioned 
earlier, giving rise to a larger volume of e~e + pairs. Moreover, no null geodesic has more than one turning point in 
its radial motion outside the event horizon. In general, a bound orbit at constant radius is considered as an orbit 
with a "perpetual" turning point. The bound unstable orbit of the photon at constant radius (the photon orbit) may 
be considered as an orbit on the verge of having a perpetual turning point. In addition, before proceeding, one way 
to described particles trajectories in the ergosphere is the following: As a particle travels with P r > or P r < to, 
say P r = 0, all local inertial frames as measured by an observer at infinity, i.e., in the BLF, are dragged into rotation 
by the KBH, irrespective of the four-momentum of the particle as it moves relative to local or neighboring inertial 
frames. (Note, particles could also be undergoing local relativistic scatterings as in the Penrose processes considered 
here.) The BLF observer measures the effects, given by P M = (P r , Pq, L, —E), of this frame dragging on the particle 



as it moves in such a strong, locally constant, gravitational environment (more on this in Section [II A). 

Now proceeding, manipulations of the results given in [5] yield, in the following form, the final four-momentum 
components of the PPP positron (e + ) and negatron (e - ) as measured by an observer in the BLF: 



E+ = e"7 c 



1/2 



(/3c.m.)rSin(9g 



COS( 



+(/3c.m.)esin0g + sin^+ + (/3 c .m.)<E> cos^+ 



7 2 
/cm. 



(/3cm.)* ( / 9 c . m .)esin0° + sin0e + 



I 7c.ro. + 1 

+ (/3 c .m.)rSin6>g + cos<//; + 



l + (/3c.„,)l 



COS0f 



+ 7c.m.C#cm.)*£+ 



1/2 



7c.ro. + 1 



(24) 



(P+)r = e w 



1/2 



7CE 



1 



(/?c.m.)r (A.m.)* COS ( 



+ [flc . m . ) e sin B c e+ sin <jf e+ 



l + (/3c.m.) 



2 7c. m. 
7cm. + 1 



Sm.0 C e+ COS <f>l+ > +7c.m.(/3c.m.)r£+ 



(25) 



(P+)e = e 



{4 ,4) 



1/2 



Tc.m. 



[ 7c. m. + 1 

+ (A:.m.)rSin0e+ COS (j) c + 



(Pc.m.)e (/3c.m.)$ COS( 



l + (A,ro.)| ' ! 



7c. m. + 1 



sin6>g+ sin0g + \ + 7 c . m .(/3cm.)e^ 



(26) 







(P+), =L + = 
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7 2 
/cm. 



-(/3c.m.)4 



[ 7c. m. + 1 

(/3 c .m.)e sin 6^+ sin^+ + (/3 c . m .) r sin 6^+ cos^+ 



l+(/3c.m.)|- 7c - m - 



7o.ro. + 1 

+7cm.C9cm.)»4 , (27) 

where, again, e", e^, e M1 , e M2 , and a; are given by Eqs. (g) through (Q), however, evaluated at the photon orbit r p h 
[Eq. (p3|)l . Note, sometimes I refer to the e~e + pairs as PPP electrons. 

The other expressions in the above four-momentum components [Eqs. ( p4| ) through (p7|)] are defined as follows 
(just as in the PCS, these definitions are derived from [5]; the results are presented here). The subscript cm. defines 
the center-of-momentum frame, and the superscript c indicates parameters measured relative to this frame. The 
cm. frame and the LNRF are related by a Lorentz transformation, with relative frame velocity defined by /3 c . m . = 
v c .m. = [(/3 c .m.)r, (Pc.m.)e, (A:.m.)*], and with the Lorentz factor of the cm. frame given by j c . m . = (1 - &l. m )~ xl2 - 
In general, 



P 7 1 + P 7 2 



(28) 



where e 7 i and p 7 i are the energy and space momentum vector, respectively, of the radially infalling photon as 
measured by a LNRF observer; e l2 and p 7 2 are these measured quantities for the orbiting photon. Using eqs. (3.90)- 
(3.92), with 8 1 i — ir/2 and </> 7 i = tt, and (2.8a)-(2.8c) of [5] in Eq. ( p8| ) above, and after some manipulation, we find 
that 



Pc.m. — [(Aj.mjrj (Pc.m.je, (/3c.m.)*J 

= {E ll +E l2 ~ojL l2 )- 1 {-E lU e"-^^/Q^,e v -^L l2 ), 



(29) 



where E^i and E l2 are the conserved initial energies of the infalling incident and target orbiting photons, respectively, 
measured in the BLF, and L l2 is the corresponding conserved azimuthal angular momentum of E l2 . The cm. frame 
energy is given by [5] 



£ 7 l£ 7 2(l — COS# 77 ) 



-E 7 i(E l2 



loL 



J2), 



(30) 



where 9 11 (= it/2") is the polar angle between p 7 i and p 7 2 in the LNRF. The polar and azimuthal angles (#g+,0g+) 
of the e + in the cm. frame, with the polar axis of the coordinate system pointed in the positive e$ direction, are 
given below in terms of the scattering angles which are found from application of the Monte Carlo method 

to the cross section for these pair production scattering events [5] : 
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arctan 
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(31) 
(32) 



where 



D l = cos 2 (9^ sin 2 6^; 
D 2 = 2 cos e c + sin d° + sin <jf + (1 ■ 
D 3 = 1 + cos 2 9% (cos 2 6% - 2 



D 3 = 1 + cos 2 6^(cos 2 
The corresponding angles for the e~ are given by 
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(34) 



jj) through ( pTj ) for the e , change the 
-> , with the new angles defined by 



Note that, to find the corresponding four-momentum components of Eqs. 
subscripts to "— ", using the fact that e c + = e'L, and let 9 c e+ — ► 

Eqs. (|33|) and ( p4| ) above. Equations ( |29| ) and ( |30| ) will be important in the discussion of the GM field in Sections |II D| 
and pl|. i.e., because of the frame dragging term u contained in these equations [cf. Eqs. through (27), (|29|) 
and (|30|)]. 

Similarly, as for the PCS photons, numerous cases of the obtained four-momentum spectra of the PPP (77 — > e~e + ) 
electrons are presented in [5] . Below, in Section III , I present only spectra that are relevant to the discussion at hand 
of how the GM field affects the Penrose produced e~e + pairs. 



D. The Gravitomagnetic (GM) Field 

In order to understand and to appreciate the origin of the asymmetry in the distributions of the scattered particles 
in the polar (O) direction, mentioned in Section |, one must look in details at the gravitomagnetic (GM) force field 
[10] — the force responsible for the asymmetry. The gravitational force of a rotating body of mass consists of two 
parts, which are analogous to an electromagnetic field. The first part is the familiar gravitational force that a body 
of rotating or nonrotating mass M produces on a test particle of mass m: 



dp 



(l_ v 2)l/2 *»' 



(35) 



where the gravitational acceleration g is produced by the gradient of the "redshift factor" e v of Eq. 
inverse e~ v is referred to as the blueshift factor [5,9]: 

g = -Vine" 

EM(r 4 - a 4 ) + 2Mr 2 a 2 A sin 2 6 , 



note, the 



2Mra 2 (r 2 + a 2 ) 



AVE 3 



cos 8 sin 6 e©. 



(36) 



The above force is analogous to the electric force field surrounding an electric charge source (i.e., like the Coulomb 
field between point charges), and for this reason is sometimes referred to as the "gravitoelectric" force. 

The second part, however, is less familiar. It is the additional gravitational force that a rotating mass produces on 
a test particle. This force, called the GM force, is produced by the gradient of /3 GM = -we$, where u> (= — <?$i/g$$ 
[7]) is the frame dragging velocity appearing in Eq. (|l]), expressed by Eq. (||). The GM force exerted on a test particle 
of space momentum p is given by [10] 



dp 

dr 



= H • p, i.e. 



GM 



dp 



(37) 



GM/ 



with 



H 



e-"V/3 GM , i.e. 



(/^GM^Ii 



(the vertical line indicates the covariant derivative in 3-dimensional absolute space). The field H is called the GM 
tensor field, and /3 GM is sometimes called the GM potential. The components of H in the LNRF are [10] 



H, 



2aMrcos9 r 
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E 2 VA 3 



a 2 )A-a 2 EAsin 2 9 



11 



2aMVA sin 9 



E 2 VA 3 
-a 2 (r — M) sin 2 9 



a 2 cos 2 6,4 - Er 



2r(r 2 +a 2 ) 



H$ r — — 



2aMVAsin6 



r z A - Er 



2r(r 2 +a 2 ) 



-a 2 (r-M) sin 2 6 \e~ v = ff e . 



(38) 



(The above definitions will become clearer in later sections.) Note that, like the magnetic Lorentz force f B = q(y x B), 
on a particle of charge q, moving with velocity v, in a magnetic field B, the GM force of Eq. ([37]) vanishes when the 
particle is at rest. 

The justification for the word "gravitomagnetic" becomes clear when we look at the GM force at a large distance 
from the event horizon (r ^S> r + ), exerted on a particle of mass m traveling with low velocity (v <C c). The GM force as 
measured by an observer at infinity [of whom I will refer to, after the introduction of Eq. (|47j), as the asymptotic rest 
BLF observer], under these conditions, is given by F GM ~ m(v x H), where H is called the GM vector field — which 
contains the same information as the antisymmetric part of the GM tensor field [10]: 

W = e* l H u = e* l e- v Vk (/3 GM ) ; ; i.e., H = e^V x & M , 



H = H e e& + H r e r 



2aM 



(' 



1 cos 2 6) sin 6e© 



2r(r 2 + a 2 ) 



cos 6e r 



(39) 



Thus, one sees that the GM force on a test particle of mass m is similar to the Lorentz magnetic force f B on a test 
particle of charge q, as given above, with the vectors H and B being analogous. In addition, the lines of force of H 
have characteristics similar to that of a dipole magnetic field. Therefore, because of resemblances such as these (more 
examples given in [10]), of the force field produced by the mass and angular momentum of the rotating black hole, to 
that of a magnetic force field, thence came the term gravitomagnetic. 

To assure a further understanding of the GM force field, one can think of it as the following. Just as the gravitational 
force causes local inertial frames to fall with acceleration g, the GM field H causes local inertial frames to rotate with 

angular velocity f2 GM = — H/2. In other words, the GM field H can be thought of as a force field that drags local 
inertial frames into rotation and, as a result, produces a "Coriolis force" m(v x H) at r S> r + . The angular rotation 
velocity r2 GM is universally induced in all bodies and neighboring inertial frames, at a given radius, i.e., independent 
of the mass m of the body, just as in the case of g. 

Moreover, Q GM is equal to the precession rate of the spin s of a gyroscope in the GM field of a rotating body [10]. 
The angular momentum of the gyroscope is changed by the torque ds/dt = (l/2)s x H, analogous to the torque jl x B 

exerted by a magnetic field B on a magnetic dipole moment jl of a "test" magnet. Also, the precession rate fi GM , 
like the acceleration of gravity g, is independent of the composition and structure of the gyroscope. It is a universal 
precession induced in all bodies, even those with vanishingly small angular velocity. 

The effect of inertial frame dragging, specifically that concerning the GM precession of gyroscopes, was discovered 
in the equations of general relativity by Thirring and Lense in 1918 [1] and, thus, is often called the Lense-Thirring 
effect. Since there has not yet been a detection of the GM field of any rotating body [10] — and to date (1997, original 
copyright date of this manuscript) this still may be the case, it would be another significant test of general relativity, 
if we find, as these calculations suggest, that the one-sidedness in the jets of AGNs is due largely to the GM field. 

I now use Eqs. (I37I ) and ( |38|) to evaluate the force experienced by a test particle under the influence of the GM 
field, as measured by a LNRF observer. In the LNRF, since the covariant components of a four-momentum vector 
{Pn = [Pr , Pe i P$ , Pt {— — £)]} are related to the contravariant components by — r/^p^, where 77^" = Minkowski 
metric components, then the space components of the GM force, exerted on a test particle, are found to be 
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(40) 



and 



(/gm) $ 



2aMrcos6 



E 2 v# 
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(41) 




E 2 \M 3 

(r 2 + a 2 ) A - a 2 EAsin 2 6 p e = H e p r - H r p e . 



(42) 



(Again, the definitions will become clearer in later sections.) Even though these GM force components are what is 
measured locally, i.e., in the LNRF, they are sufficient for our purposes to describe what we observe at infinity, i.e., 
in the BLF, since the space momentum components in the LNRF (p r ,po,p$) and the BLF (P r ,Pe,P$) are related 
in a linear fashion (see the transformations between these frames given by eqs. (2.7) and (2.8) of [5]). Notice that 
the radial component (f GM ) r and the polar component (/ GM )e [Eqs. ( fi"o"| ) and (f4lj)l, exerted by the GM field, are 
proportional to the azimuthal momentum of the particle. This means that the absolute magnitudes of these force 
components will increase as the energy, E', of the scattered particle increases, since (= L') and E' increase linearly 
[5]. The azimuthal component (/ GM )$ of Eq. ( pf2] ) will act according to the radial and polar momenta of the incident 
and scattered particles, exerting a force in the azimuthal ($) direction. This GM force component could be important 
in the polar collimation effects produced by the rotating black hole [25,26]; however, it is beyond the scope of this 
manuscript to discuss it in that context here. 

Displayed in Figs. 2(a) and 2(b) are the ratios (f GM ) r /p<s> and (/ GM )e/p*, respectively, as functions of the polar 
angle O, which will be compared to the polar angle of escape for a Penrose scattered particle. The angle 9 varies 
from to 7r, with the equatorial plane located at ir/2. For a given positive as opposed to negative p$ for 
retrograde motion, notice that (/ GM )r of Fig. 2(a) reaches a maximum in the equatorial plane, and increases with 
increasing r away from the event horizon, within the range of r considered. In Fig. 2(b), we notice that (/ GM )e is 
antisymmetric, above and below the equatorial plane, and reaches negative and positive maxima at O ~ 60° and 
~ 120°, respectively. Moreover, the maximum value of the GM force component (/ GM )e is larger than that of 
(/ GM ) r by factors of ~ 10.5 and 16.4 at r m b and r p h, respectively, and the strength of (/ GM )e decreases little with 
increasing r [cf. Fig. 2(a) and 2(b)]. Figure 2 will be discussed further in the following section. 



III. DISCUSSION: THE GM FIELD AND THE PENROSE PROCESSES 



In this section, I use the results of the four- momenta of the Penrose scattered particles, found in [5], and presented 
again in Eqs. (O) through JTs| ) and j24| ) through (|27|) showing specific details, to see what internal role the GM force 
of Eqs. (|40| ) through fl42j) has in determining the spacial distributions of the Penrose scattered escaping particles. In 
the following, one must bear in mind that the GM force has already been accounted for in these Penrose processes 
(as done in the model calculations of [5]) through the Kerr metric tensor of Eq. (|l|), whose function is to describe 
the geodesies of the particles, on which the four-momentum is defined. That is, just as we do not have to calculate 
separately the gravitational force on the scattered particles [Eq. fl3q)], we do not have to calculate separately the GM 
force on the scattered particles [Eq. (p7|)], because both of these gravitational forces are inherently tied to the overall 
scattering process through the Kerr metric. The GM force, however, is given here separately only to show how it 
relates to the final four-momentum vectors of the scattered particles. Therefore, when looking at this force separately, 
in some instances, the precise relationship of the GM force with a distribution of scattered particles may not be clear. 
With this in mind we now proceed. 
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A. The Gravitomagnetic Field and PCS 



Beginning with the PCS, displayed in Figs. 3—7 are scatter plots of the momentum components {Pp h ) r and (P ' h )e 
[see Eqs. ( |l6|) and (0)], and the corresponding polar angles 6 ph [to be compared with in Eqs. ( pZ)|) through (42)] of 
the escaping PCS photons. The polar angle (defined here as the angle of escape [3,4]) is given generally for a particle 
of type i by [25] 



=F arccos 



-T+ VT 2 -4S77 
2S 
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where 



S = a 2 (£-E 2 ), 
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the negative and positive signs are for particles traveling above or below the equatorial plane, respectively. The polar 
angle of Eq. (|4^) is derived from Eq. ( |l2"| ) by letting Pq — > 0. (Note, such scatter plots as displayed in the figures 
presented here result from the application of the Monte Carlo method, which yields individual scattering angles and 
four-momentum components for the scattered particles.) 

Upon comparing Figs. 2(a), 3, and 6, we observe the following. For PCS at r = r m b, between radially infalling 
photons (-Eph < 1 MeV) and tangentially equatorial orbiting target electrons [Figs. 3(a)— 3(e)], the spectrum appears 
symmetric in (Pp h ) r for the lowest incident photon energy used: E p h — 0.511 keV. What role the GM field plays in 
giving rise to this symmetric distribution is unclear [cf. Fig. 3(a)], since such a distribution is expected in general. 
In Figs. 3(b)— 3(e), for the higher energy incident photons, it appears that GM force component (f GM ) r of Fig. 2(a) 
has little effect on the distribution of escaping photons. It seems that the distribution of the scattered photons in 
the radial direction is, in these cases, dominated by the momentum forces involved in the scattering process itself. 
That is, as £" p h is allowed to increase, while E e is held constant at 0.539 MeV, the preferred direction of the scattered 
photon is radially inward — such behavior is to be expected. In addition, even though some of photons are scattered 
in the positive radial direction, the scattering behavior is inconsistent with what would be expected if (f GU ) r were 
important, namely, one would expect for the particles with ph ~ tt/2, and large (-P ph )$, and thus large E' ph , to have 
positive radial momenta [cf. Figs. 2(a), 3(b)— 3(e), 6(b)— 6(e), and Eq. flic])]. 

In the case of PCS by nonequatorially confined target electrons [Figs. 3(f)— 3(i), with 6.145 M 2 m 2 < Q e < 
614.5 M 2 m 2 ] — i.e., as the energy- momentum of the target particle is allowed to increase, we find that the num- 
ber of photons with outward radial momenta increases with increasing E' ph [implying increasing (-Pp h )$]- When one 
compares the higher energy photons and their polar angles of escape, displayed in Figs. 6(f)— 6(i), with Fig. 2(a) and 
eq. (40), it appears that (/ GM )r may have some importance in the behavior causing the number of photons scattered 
with positive radial momenta to increase [cf. Figs. 3(f)— 3(i)]. Such behavior could be due, in part, to the fact that 
(f GM ) r is large on the particles with large (-P ph )$ (i.e., large E' ph ) and with 6 ph near ~ 7r/2; this is consistent with 
what we find in the cases of Figs. 3(f)— 3(i) and corresponding angles 6(f)— 6(i). The other part responsible for the 
behavior could be that, as E e become 3> -Eph, with increasing Q e , backward scattering [(PpiJr > 0] becomes more 
effective: cf. 3(d) to 3(f) and 3(g), cf. 3(c) to 3(h) and 3(i), showing progression towards symmetry as E e becomes 3> 
E p h- We will return to this discussion after the introduction of Eq. (|47|). 

Before proceeding, note the following: the nonequatorially confined (spherical-like orbiting) target electrons dis- 
cussed here, and target photons discussed in the next section, for a given distribution, can have either Pq < or 
Pq > 0, as the particle passes through the equatorial plane going above or below, respectively. In the figures shown 
here, one-half of the distribution of 2000 targets are given the negative Pq value and other half are given the cor- 
responding positive Pq value (this was not done in [5]). Such a distribution yields what would likely be seen by an 
observer at infinity. Also, note that, consideration has been given only to collisions between infalling incident particles 
and bound orbiting target particles, thus, assuming that collisions between neighboring target particles are negligible. 

Continuing, now the discussion centers on the various distributions of PCS photons as they escape in the positive 
and negative e© directions. Upon comparing Figs. 2(b) and 4—7, we observe that the polar angles 9 ph of the 
PCS photons, after scattering off equatorial target electrons [Figs. 6(a)— 6(e)] and nonequatorially confined target 
electrons [Figs. 6(f)— 6(i) and 7(a)— 7(d)], are consistent with the GM force component (/ GM )e [Eq. ([ll])] acting on 
the scattered particles. The effects can be seen in the polar coordinate momenta (PLyJe, displayed in Figs. 4 and 5 — 
corresponding to cases of Figs. 6 and 7, respectively. Looking closely at the distributions resulting from the equatorial 
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target electrons, the PCS photons go from the general expected symmetric distribution, for the lowest monochromatic 
incident photon energy used: E p h = 5.11 x 1CP 4 MeV, to being strongly asymmetrical or one-sided, for the highest 
energy used: E p h = 1 MeV, with E e ~ 0.539 MeV being the same for all the cases; cf. Figs. 4(a)— 4(e). Overall, we 
find that, as the energy of the incident photon E p h is increased, more and more of the photons are scattered with 
{P ph ) r < 0, (P ph )e > [i.e., with 9 ph > 90°; cf. Figs. 3(a)-3(e), 4(a)-4(e), and 6(a)-6(e)]. Apparent reasons for 
such distributions are explained in the following. 



Recall that the magnitude of the GM potential is | /3 GM |= lo (see Section [ID ). By examining the angle of Eq. fl2l|), 
we can see how this GM potential, producing the GM tensor field [see Eqs. (37|) and (38)], causing the frame dragging, 
may alter the incoming angle 8 R h of the incident photon, and subsequently the outgoing angles, 9 R h and (/>p\, of the 
scattered photon, given by Eq. (E3) and Eq. (3.39) of [5]: 



b R ' 6 R - arccos [ cos ^ ~ cos <h cos <h 
Ph ~ 0ph C { sin «ft sin 0* 



(44) 



respectively, where 4>p\, the important angle here, gives the direction of the scattered photon, above and below the 
equatorial plane. [Note, the angles 9 R h and <fi R h define the space momentum components as described by a unit 
four- vector tangent to the coordinate lines at the event (r, 0,$,t), in the ERF [5]; this particular LNRF coordinate 
system is centered on the electron with the pole in the direction of positive e^.] So, in essence, what we are finding is 
that, the GM force close to event horizon, whose presence is seen by the appearance of the frame dragging velocity u 
[cf. Eqs. (|2^) and (0)], exerts a force on the moving particles, causing them to behave in various general relativistic 
ways (uncharacteristic of Newtonian physics), as discussed further below. 



Specifically, for the cases of Figs. 4(a)— 4(e), and corresponding Figs. 6(a)— 6(e), Eq. (21) shows that 9 ph will not 
change for PCS off monochromatic distributions of 2000 target electrons, at r m b, if E e (~ 0.539 MeV) is held constant 
for the different cases while E^ is allowed to increase, increasing from 5.11 x 10~ 4 MeV to 1 MeV. But, increasing 
i? p h will, however, change the scattering angle 8 R , subsequently, still, altering the scattering angles 9 R h and 4> p h [cf- 

Eqs. ( p2f ) and (p4|)]. So, how does S R relate to the frame dragging velocity uil The angle S R is determined, as stated 
earlier, from the Klein-Nishina cross section: the effective area of an electron producing a scattering event in which 
a photon is emitted at particular polar and azimuthal angles, S R and a R , respectively, of a spacetime coordinate 
system centered on the electron, with the pole in the direction of — e r (i.e., of the initial photon direction). The angle 
S R depends on the Lorentz factor 7 e of the target electron and the initial incident photon energy I? p h [see Eq. (p0|)]. 
Equation ( |l9|) reveals the direct relationship between 7 e and oj, and, thus, 5 R . 

Investigating further, the angle 4>p\, through the space momentum vector of the scattered photon in the ERF [5]: 



(45) 



(PpL)r 


= Pph sin' 


?ph cos 


& 


(Pph)e 


= P% sin ' 


?ph sin t 


bph 




= Pph cos 


n R' 
Cphi 





gives, in general, in the four quadrants relative to the equatorial plane, that 

I : < <j) R h < 90°, II : 90° < <j) R h < 180°, 

(P%)r > 0, (p£) P < 0, 

bpt)e > 0, (pg) e > 0, 

90°< 6p h < 180°; 90° < 8 ph < 180°; 

III : 180° < 4> R i < 270°, IV : 270° < < 360° 

(P%)r < 0, (pg) r > 0, 

«h)e < 0, (pg) e < 0, 



0< 6 ph < 90°; < 6 ph < 90°. (46) 



We will use the quadrants of Eq.(|46|), in the discussion below, in an effort to further understand the distribution of 
PCS photons. (Note, the Lorentz spacetime transformations for eq.(45), from the ERF to a general LNRF observer, 
are given in [5] [eq. (3.46)]). 
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Upon comparing Eq. ( |46| ) and Figs. 4 and 6, we find the following. In Fig. 4(a) [see also Fig. 3(a)], the photons are 
scattered symmetrically into the four quadrants. Increasing E p h causes the highest energy photons to be scattered 
below the equatorial plane, and the appearance of a void in the lowest energy regime; cf. Figs. 4(a)— 4(e). These 
figures show how asymmetry in the preferred positive e© direction, of (P p h)e>, increases as E p ^ is increased. The 
highest energy photons are progressively scattered into quadrant II, consistent with (/ GM )e [see Eq. (41) and Fig. 2(b)] 
acting below the equatorial plane on these scattered particles, proportional to increasing (P ph )<s> (i.e., increasing E' ph ), 
with 90° ^ 6 ph ^ 180°. On the other hand, the lowest energy photons are progressively scattered into quadrant III, 
consistent with (/ GM )e acting above the equatorial plane, producing a smaller GM force on these particles, and thus 
a smaller magnitude for (P ph )e, with ^ 9 ph ^ 90° [cf. Figs. 6(a)— 6(e)]. The most probable explanation for these 
distributions as related directly to the GM field, although not apparent from Eqs. ( f40| ) through ([l|]) — which give the 
GM force as measured by a LNRF observer, will be given below. 

So, overall, the equatorial target electrons scatter the photons preferentially in the direction of positive e@, as the 
energy-momentum of the incident photon is increased. The asymmetry seen in the scattered particle distributions 
described above, at these relatively low incident and target particle energies, while E e (~ 0.539 MeV) is held constant, 
can be viewed as the "signature" of the GM field acting on the infalling [(P p h)r < 0] incident photons, and the scattered 
photons, while increasing E p h (i.e., increasing |(-P p h)r|)- The resulting effects can be seen in the momentum component 
of Eq. ([17j), plotted in Figs. 4(a)— 4(e), in which the angles contained therein are altered by certain combinations of uj 
and the energy-momenta of the incident and scattered particle energies. The ratio of the number of escaping photons 
with {P ph )e > to those with {P ph )e < (= e*, for particle of type i) range from e ph ~ 1.07, ~ 1.46, ~ 2.42, ~ 3.42, 
up to ~ 4.99, for Figs. 4(a)— 4(e), respectively. Thus, as stated earlier, it appears that reflection symmetry of the 
Kerr metric, above and below the equatorial plane, is broken due to the GM force field. 

Note, in general, since is positive for all of the escaping Penrose scattered particles, the antisymmetric 

component (/ GM )e [see Eq. ( fil| ) and Fig. 2(b)] will exert a force in the negative e@ direction for O- < 90° and in the 
positive e© direction for > 90°; such could aid in the acceleration of the jets of AGNs. Namely, this GM force 
component has the effect of blueshifting the escaping photons to higher energies, while increasing the Lorentz factor 
of the escaping electrons. 

Now we will consider the PCS by nonequatorially confined target electrons, with the initial monochromatic photon 
energy held constant at -Eph = 5.11 x 10~ 4 MeV, while E e is allowed to increase from the equatorially confined value 
E e ~ 0.539 MeV, Q e = [Fig. 4(a)] up to E e ~ 5.927 MeV, Q e = lM.5M 2 m 2 e , for the different cases [cf. (P ph )e and 
corresponding angles of escape O ph of Figs. 5 and 7, respectively]. The asymmetry proposed to be due to the GM force 
(indicated by the prominent void) is quite apparent in the different cases [cf. Figs. 5(a)— 5(c)], i.e., before, it seems, 
momentum forces inherent to the scattering processes, because of increasing Q e , and the final effects of the GM force 
acting on the scattered particles, dominate over the initial effects of the GM force, thus, causing symmetry to reappear, 
as seen in Fig. 5(d). For the nonequatorially confined target electrons, I did not derive expressions as done for the 
equatorial targets [Eqs. (|l5| ) through (p2|)], showing the direct relationship between the momentum components and 
the angles producing the asymmetry. Nevertheless, it appears that the asymmetry is caused, in general, by the same 
mechanism, i.e., the GM force alters the incoming and outgoing angles of Eqs. (|2^) and (H), through the presence 
of ui. To further compare, the ratio of the number of escaping photons with (P ph )e > to those with (P ph )e < 
range from e ph ~ 1.07 for Fig. 4(a): Q e = to e ph ~ 1.54, ~ 2.13, ~ 1.29, ~ 1.02, for Figs. 5 (a) -5(d): Q e ^ 0, 
respectively. The specific cause of the behavior for these scattered photon distributions, varying from symmetric to 
asymmetric [scattering preferentially into quadrants I and II; see Eq. (46)], then back to symmetric (scattering equally 
into the four quadrants), is not readily accessible from the GM force components [Eqs. ( [lo| ) through (0)], but, again, 
the presence of the frame velocity in the scattering angles [Eqs. (22) and ( |44| )] assures us that the GM force field 
has a role in determining the trajectories of the photons. We will return to these distributions below, where we will 
find that the asymmetry seen in the scattered particle distributions, at the relatively low incident and target particle 
energies [Figs. 5(a)— 5(c)], while E p h (= 5.11 x 10 -4 MeV) is held constant, can be viewed (as in the case of the 
equatorially confined target electrons described above) as the signature of the GM field acting on the initial infalling 
photons, and the scattered photons. [Note that, the highest energy photons are scattered with polar escape angles ph 
approximately equal to the polar angle indicating the maximum or minimum latitudinal angle of the spherical-like 
orbiting target electrons (~ 29.5° relative to the equatorial plane), given by Eq. (|43|); cf. Fig. 7(d).] 

Deviating briefly with a discussion to help us to understand better the above scattered particle distributions and 
those to follow: we do not yet have analytic expressions for the force components [like Eqs. (^) through (Q)] 
describing the behavior of the GM field on a test particle near the event horizon as measured by a local BLF (LBLF) 
observer. The reason for this is because of the complexity in deriving an analytic expression for the GM tensor field 
H inside the ergosphere in the LBLF, i.e., due to the severe dragging of local inertial frames, in which the time 
coordinate basis vector (d/dt) of the Kerr metric changes from timelike to spacelike [7]. Algebraically, this complexity 
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arises in the partial derivatives of the Kerr metric components, contained in the nonzero Christoffel symbols 1^ 
(sometimes referred to as the affine connections), when evaluating the covariant derivatives in the general expression 
for Hij [cf. Eq.(^)]. Therefore, until an exact expression for H can be found, we can only approximate its physical 
nature based on the information we have at hand, and what these Penrose scattering processes are telling us, about 
the behavior of the GM field close to the event horizon. Now, we know that the GM field at radii 3> than the event 
horizon (r — r + ) appears as a dipolar- like field as measured by a BLF observer. This dipolar- like field will however 
be distorted near the event horizon as measured by this distant observer, due to large velocities and strong gravity, 
but not fully destroyed [10]. It is possible that the preference of the particle being scattered into the positive e© 
direction is a display of this distorted dipolar-like field. The most probable explanation for the preference in the 
positive e© direction is that, near the event horizon in the BLF, the GM force field lines are frame dragged into 
the direction of the rotating KBH, i.e., positive azimuthal (e#) direction, which results in H, the GM tensor field, 
acquiring additional cross terms: H t q and H® r [cf. Eq. (|38|)]. This is supported by the acquired nonzero T^ 's in the 



y 

covariant derivatives — appearing in the component notation for H [see Eq. (B7|)], which in general gives rise to nonzero 
H t q and ffer components. In this case, the GM vector field, say H, will acquire an azimuthal component ff*, such 
that H = H r e r + ff e e© + i?*e* [cf. Eq. (p9|)]. Then, by the analogy of the approximation in the asymptotic rest 
(or nonrotating) BLF: F GM ~ P x H, for r 3> r+, and using Eq. we find that the GM force on a test particle as 
measured by a LBLF (inertial framed dragged) observer, for r ~ r+, can be expressed by the following: 

F GM - H • p 

= [( F GM.)r> C^gm)©) C^gm)*] 

= {H re p e + ff r $p*)e r + (H Br p r + ff e $p*)e© + (ff*rP r + ff$ep e )e* 

oc (ff*P e - ff e P$)e r + (iTP* - H*P r )e@ + (ff e P r - ff r P e )e*. (47) 

Note, the oc indicates that we are assuming (1) the space metric tensor g 3 , multiplying the covariant momentum 
components has been defined by a local orthonormal (Lorentz) tetrad h£ so that, as usual, g^ v — h/ffiprf*" with 
rf-P being the Minkowski metric components; and (2) the space momentum transformations from the LBLF to the 
global BLF are standard [7], i.e., "somewhat" similar to transformations from the LNRF to the global BLF [5]. It 
is conceived that the above assumptions are at least valid for our analyzing purposes. Again, I remind the reader 
that the GM force presented in Eq. (Et|) is derived here only to show how the particles might be affected by the GM 
field, but these effects, shown in the resulting distributions of the escaping particles as measured by a BLF observer, 
are intrinsically incorporated into the calculations through physical processes occurring in the Kerr metric spacetime 
geometry. The components of the vectors H and H in Eq. ([47]) are defined by the components of the GM tensor: 



ff r 


= ffe$, 


H r = 


— ff$e 


H e 


= ff$ r , 


ff e = 


— ff r $ 


ff* 


= H r Q] 


ff* = 


— Hq t 



(48) 

[cf. Eqs. (ff), ©, and © through ©]. If the GM tensor is antisymmetric [sec Eq. ©], then H r = H r , H e = H , 
H* = H* 

From a geometrical analysis of the characteristics of the antisymmetric GM tensor components H in the asymptotic 
rest BLF (ABLF) [see Eq. (|39"|)] and the components of H in the LNRF in which the frame dragging is "taken out" 
[see Eq. (pq)], and the probable condition in the LBLF, that the GM tensor field lines are frame dragged in the 
azimuthal direction, we can deduce that the GM tensor in the LBLF will be somewhat similar to that in the LNRF, 
in respect to the dipolar-like characteristics, with differences, in general, depending on the frame dragging velocity u> 
[cf. Eq. (p7|)]. The tensor components of Eq. (p8f), when evaluated at r ~ r + , reveal to us the following: 



H r < and H r < for 9 < 90°; 
H r > and H r > for 9 > 90°; 
H e > for 9 < 90° and 9 > 90° 
H e < for 9 < 90° and 9 > 90° 
ff* < for 9 < 90° and 9 > 90° 
ff* > for 9 < 90° and 9 > 90° 



(49) 



cf. Eqs. ( [47| ) and d48|). The additional components included in Eq. (|49|), ff* and ff*, are found from deductive 
reasoning as follows. We know from a general direct evaluation of the covariant derivative of /3 GM that H t q = Hq t . 
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Therefore assuming that H r Q = Hq t > (as would be expected because of the frame dragging), then by definition 
[Eq. (fl8l)l and the results of the Penrose scattered particle distributions, we can deduce the general signs of if* and 
if* in, above, and below the equatorial plane. We will use the relations of Eq. ([l9]), as we continue in the discussion 
below. 

Note, the differences between the frames of reference considered here are the following. An observer in the BLF 
measures coordinates r, O, $, and t, relative to the KBH. The BLF is an inertial frame (i.e., it is a freely falling 
coordinate system) . An observer in this frame measures spacetime events from two points of view: the asymptotic rest 
(or flat) inertial frame (events occurring at r ^ r + ) — which I refer to as the ABLF, and the global inertial frame (events 
occurring at r ~ r + ) — which I refer to as the BLF or global BLF. The measured four-momenta [Pp = (P r , Pq, L, —E)], 
describing the trajectories of the particles [cf. Eqs. J14 ) through (Q)] , are the same parameters for both these observers, 
although the immediate surroundings of the particles are very different. The BLF observer at r ~ r + measures the 
effects of frame dragging, whereas the ABLF observer does not, yet both observers are located at infinity. For example, 
the GM tensor field as measured by the ABLF observer will have prominent dipolar- like features [cf. fl39|)], features 
that will be distorted as measured by the LBLF and BLF observers. An observer in the LBLF "sees" this distorted 
field, whose force lines are frame dragged into the azimuthal direction, and measures its effects on moving test particles 
[cf. Eq. (p7|)], but the ABLF observer cannot measure these near event horizon effects; it is only through the "eyes" 
of the LBLF and BLF observers, as evident in these Penrose process distributions, that the observer at infinity can 
measure these ergospheric effects. Now, on the other hand, the LNRF observer is a local inertial frame rotating with 
the geometry [7], experiencing in a sense no frame dragging. This spacetime coordinate system is at rest relative to 
the frame dragging, and revolves with velocity u> relative to the observer at infinity. Physical processes described in 
this local Lorentz frame are "simpler" because the frame dragging is canceled as much as possible. The observer in 
the LNRF measures the force that the local gravitationally distorted dipolar-like GM tensor field exerts on a moving 
test particle, however, free of frame dragging effects. Therefore, the LNRF has its limitations in explaining the local 
effects of the frame dragged GM force field; for this reason, the LBLF is needed. 

Equation ( [l7j) shows that the GM force on a test particle in the LBLF and the BLF (wherein we are assuming the 
most probable condition) will acquire additional terms in the radial and polar components, different from that in the 
ABLF and the LNRF. [Note, in the LNRF, Eqs. (^) through (Q) have the same general form as Eq. J47), however, 
with fi* = if* = 0.] These additional terms involving if* and if*, as we shall see, can explain the asymmetries that 
we observe in the Penrose scattered particle distributions presented here. For a preliminary example, if P r is less than 
zero, and allowed to increase in magnitude for the different cases considered, as for the initial radially infalling photon 
[(Pph)r < 0, (Pph)e = (P P h)* = 0] displayed in Figs. 3(a)— 3(e); 4(a)— 4(e); 6(a)— 6(e), with the target particle having 
relatively low energy-momentum, Eq. ( B^ ) shows that (P GM )e will exert a dominant force on the photon, increasing 
in the positive e© direction, and thus, initializing the asymmetry eventually seen in the high energy PCS photons 
[cf. Figs. 4(a)— 4(e) and 6(a)— 6(e)]. Now, with this initialized asymmetry, the final distributions are consistent with 
(P GM )e, again as given by Eq. (J47|), acting directly on the PCS photons: In addition to this force component being 
positive for photons scattered with (Pp h ) r < and H r > [i.e., Op h > 90°; cf. Eq. (|49|)1, this force component will 
be positive for photons scattered with {P^ h ) r > and H r > [i.e., S' ph > 90°] provided ff r (P^ h )$ > F*(P^ h ) r , 
as for the photons with large L' ph or small (Pp h ) r — otherwise this force component will be negative. For (Pp h ) r < 
and H r < (i.e., 6p h < 90°), (P GM )e will be > for sufficiently small L' ph or sufficiently large |(P ph )r| — otherwise 
(P GM )e will be < 0. This behavior is consistent with Figs. 4(a)— 4(e). On the other hand, if the energy-momentum of 
the infalling incident photon is relatively small compared to that of the orbiting target electron [in which, say (P e )$ is 
allowed to increase for different cases, with (P e ) r = 0], Eq. ( fi7j ) shows that (P GM )e: exerted on the scattered photon, 
will dominate, creating or recreating symmetry in the polar direction. We find this behavior to be consistent with 
what we observe in the PCS photons of Figs. 4(a), 5, 6(a), 7. More on this later. 

So, in general, Eq. ([!?]) shows that it is the if* component of the GM tensor field in the LBLF, acting on radial 
component of the infalling incident particles, that subsequently causes the polar angle 0[, of the scattered particles of 
type i, to be > 90° (i.e., Pq > 0). Once the preference of the scattering angle OJ > 90° has been established by the GM 
force as measured by a LBLF observer, we find that the scattered particle trajectories are somewhat consistent with 
the GM force as measured by a LNRF observer [ Eqs. ( f4(i| ) through (|4^)] — however, the tra ject ories cannot be fully 
explained in this frame. Therefore, in the remainder of this discussion, we will now use Eq. (|l7|), the proposed, most 
probable GM force measured by a LBLF observer, which includes the use of Eqs. (^) through ([l2|) as measured by a 
LNRF observer, to further show the effects that the GM field has on these Penrose scattering processes — remembering 
that the BLF observer, i.e., the observer located at infinity measures both near event horizon events (r ~ r + = BLF 
or global BLF) and faraway events (r » r + e ABLF). 

Before proceeding, we return to the cases of Fig. 5 to shed some additional light on these distributions using Eq. (|47|). 
The radial distributions of these PCS photons remain nearly symmetrical for the different cases (not displayed here), 
looking very much like the case of Fig. 3(a), where there is a slight preference for the number of particles escaping 
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with (P^ h ) r < 0. Using Eq. @, Eq. @ shows in general that for H r > (i.e., 9 4 > 90°) and P r < 0, (P GM )e 
is > 0. But for P r > and H r > (i.e., 9^ > 90°), (P GM )e is > or < for sufficiently large or small azimuthal 
momentum Pj> (= L), respectively. On the other hand, Eq. ( |47j ) shows in general that for H r < (i.e., 9^ < 90°) 
and P r > 0, (P GM )e is < 0. But for P r < and H r < (i.e., 9, < 90°), (P GM )e is < or > for sufficiently 
large or small azimuthal momentum P$, respectively. Thus, the above analysis shows possible symmetry in the polar 
direction if there is symmetry in P' r of the scattered particles distribution. However, if there is asymmetry in P r , say, 
the particles are scattered preferentially with P' r < 0, there will be asymmetry in the polar direction: For P' r < 0, the 
component (P GM )e acting on the particles with 9- > 90° will be > 0, but (P GM )e acting on particles with 9- < 90° 
will be > for particles with \H r P^\ < |.ff*P^| (implying small L', low E' or large \PJ.\) and < for particles with 
\H r P^\ > |P*P^| (implying large L', high E' or small \P^\). This behavior of Eq. j47| ) is sufficient to give a general 
explanation for the polar asymmetries occurring in these Penrose scattering processes if there is asymmetry in the 
scattered particle spectral distribution of P' r [cf. Figs. 3(a)— 3(e) and 4(a)— 4(e)]. Yet, if there is symmetry in P' r [or 
near symmetry — as we find for the cases of Figs. 5], Eq. fl47| ) shows that (P GM )e acting on the scattered particles 
should produce symmetry in the polar direction. However, because of the initial asymmetry, favoring the positive e© 
direction, caused by the GM force component (P GM )e acting on the radial component (P p h) r of the infalling incident 
photons [cf. Eq. (47)], the scattered particle distributions, at sufficiently low orbital energy- momentum for the target 
particle, will not reflect the expected mirror symmetry, above and below the equatorial plane — contrary to what one 
would expect in Newtonian physics or in the weak gravitational regime (r 3> r+); this can be seen in Figs. 5(a)— 5(c). 
That is, according to the first term in (P GM )e [cf- Eq. (f|7|)] and the description above for small V (or for the satisfying 
of |PPP4| < -PrDi this initial asymmetry in the +e© direction is amplified when (P GM )e acts on the subsequent 
PCS photons. But, as Q e of the nonequatorially confined target electrons is allowed to increase for the different cases, 
i.e., increasing the magnitudes of (P e )e and L e , and subsequently (P ph )$, the first term in (P GM )e dominates over 
this initial asymmetry, causing the reappearance of the symmetry we observed in Fig. 5(d), thus, showing us the final 
result of (P GM )e acting on the scattered high energy photons, increasing in magnitude with increasing L' ph . So, in 
essence, (P GM )e gives the escaping particles with large azimuthal component momentum a boost into opposite polar 
directions, and is greater in the +e© direction for particles with (P ph ) r < [cf. Figs. 3(a)— 3(e) and 4(a)— 4(e)]. 

Now proceeding with the discussion after the brief deviation, the momenta (P ph )e 01 Figures 4(f)— 4(i), from PCS of 
2000 radially infalling photons (for each case) by nonequatorially confined target electrons, show presence of the GM 
force acting on the particles. Effects of (P GM )e [of Eq. fl47|)] acting on (P p h)r of the initial photons, P p h = 0.15 MeV 
[Figs. 4(f), 4(g)] and P p h = 0.03 MeV [Figs. 4(h), 4(i)], can clearly be seen, by the presence of the most prominent 
regions void of particles [cf. Figs. 4(d) and 4(e)]. This feature (viewed as the signature of the GM field acting on the 
infalling incident photons, and the scattered photons) is somewhat permanent at these initial energies, in the sense 
that (P GM )e acting on the PCS photons resulting from the target nonequatorially confined electrons, of initial energy 
E e ~ 1.297 MeV, E e ~ 5.927 MeV, and E e ~ 11.79 MeV [Figs. 4(f)— 4(i)], the "symmetry producing agent," will 
not dominate over the effects producing this feature [as in Fig. 5(d))]. Once the preference for 9 ph > 90° has been 
established by (P GM )e acting on the infalling incident photons, as described above, the additional asymmetry seen 
in Figs. 4(f)— 4(i) is consistent with (P GM )e acting on the escaping photons, which are scattered predominantly with 
(P ph ) r < [see Figs. 3(f)— 3(i)]. According to Eq. @, for (P ph )r < 0, (P GM )e will be > (or in the +e© direction) 

for all the photons scattered with 9 ph > 90°, but those with 9 ph < 90°, (P GM )e will be > if \H r L' ph \ < |P*(P ph )r|, 
as for the lowest energy PCS photons; otherwise (P GM )e will be < 0, i.e., for photons scattered with sufficiently larger 
energies. This is the behavior we observe in Figs. 4(f)— 4(i) [cf. also corresponding escape angles of Figs. 6(f)— 6(i), 
respectively]. For example, notice how (P GM )e acting on the PCS photons resulting from the nonequatorially confined 
target electrons (the symmetry producing agent) of increasing orbital energy-momentum, for the different cases, causes 
more and more high energy PCS photons to be emitted in the negative e© direction, and, thus, eventually causes 
symmetry to somewhat reappear except for the GM signature [cf. Figs. 4(h) and 4(i) to 5(d)]. Moreover, according 
to Eq. @, as E e increases and more and more of the photons are scattered with (P ph )r > 0, for 9 ph < 90°, (P GM )e 
will be < 0, assisting in reestablishing the symmetry [cf. Figs. 3(f)— 3(i)]. Note, the ratio of the number of escaping 
photons with (P ph )e > to those with (P ph )e < 0, e ph , range from ~ 2.41 for Fig. 4(f) to ~ 1.34, ~ 1.16, and ~ 1.08 
for Figs. 4(g)— 4(i), respectively. 

Astrophysically, these results mean that an observer at infinity will observe "photon jets," collimated concentric 
the polar axis, ranging from symmetric to asymmetric, depending on the accretion disk properties: which populate 
the target particle orbits and supply the infalling incident photons. So, overall, in summary, it appears that the angle 
</> ph , measured perpendicular to the equatorial plane, in the ERF (a particular LNRF observer), given by Eq. (Q), 
is altered by the presence of uj, the frame dragging velocity, such that certain combinations of increasing P p h and/or 
E e (i.e., increasing the energy- momenta of the incident and target particles) cause the scattered particle distributions 
to vary from symmetric to asymmetric, with the preference in the direction of positive e©, then back to very nearly 



19 



symmetric. The behavior of these distributions is in some degree consistent with the GM force as measured by a 
LNRF observer [Eqs. ( fiof ) through (^) ], but this behavior cannot be fully explained in this frame. The most probable 
explanation, that near the event horizon the GM held lines are frame dragged in the direction as measured by 
a local frame dragged inertial observer [leading to the derivation of the general expression Eq. (^) for the GM force 
as measured locally (at r ~ r + ) by a BLF observer], can indeed explain this behavior. We shall hnd this to be true 
also for the PPP discussed below. 



B. The Gravitomagnetic Field and PPP (77 — > e e+) 

Now considering the PPP (77 — > e~e + ), displayed in Figs. 8 and 9 are scatter plots of (P=p)e [Eq. (p6|)] and the 
corresponding polar angles 8 T [Eq. flls])], respectively, of the escaping e~e + pairs for various cases, after 2000 PPP 
events at the photon orbit for each case. In general, the force component {F gm )q of Eq. (f47|), using Eq. ([l9"|), reveals 
to us that infalling particles (P r < 0) and outgoing particles (P r > 0) moving along the equatorial plane (@i = 90°) 
will experience a force (F GM )e > and (F gm )q < 0, respectively. The GM force component {F gm )q > acting on the 
infalling incident photon of energy E 7 i = 0.03 MeV is expected be somewhat significant in initializing the asymmetry, 
above and below the equatorial plane, favoring the positive e© direction (6^ > 90°), eventually seen in the scattered 
particle distributions, in the relatively low energy regime, as indicated in Fig. 4(c). However, this asymmetry is not 
apparent in the case of Figs. 8(a) and 9(a). In this case it appears that, this initialized asymmetry is balanced (or 
canceled) by (-F GM )e acting on the outgoing e~e + pairs — at these relatively low incident and target particle energies 
(£ 72 ~ 13.54 MeV), for 6 T < 90°, and if the condition is satisfied that \H r L T \ < \H*(P^) r \ for 9 T > 90°, (F GM ) 
will be in the negative e© direction: Such could very well balance or cancel out the effects of (F GU )& > acting on 
the initial infalling incident photons. Note that, the only PPP electrons satisfying the escape conditions are those 
with (Pzp)r > 0; electrons with (P=p) r < fall into the black hole since the PPP (77 — ► e~e + ) takes place at the 
photon orbit [5]. Now, on the other hand, at larger energies, as in the cases of Figs. 8(b)— 8(d) and 9(b)— 9(d), the 
initialized asymmetry, produced by (-F GM )e acting on (P 7 i)r < (of the infalling incident photon), becomes less and 
less balanced as L T increases (implying increasing E T ), as the condition \H r L T \ > \H (P^) r \ becomes satisfied, at 
which (F' GM )e behaves like an antisymmetric function about the equatorial plane, with (-F GM )e < for < 90° 
and (F GM )e > for 6 T > 90°, for (P^) r > 0, superimposed with the established initial preference for 9 T > 90°. So, 
underlining the cause of the symmetry and asymmetry we observe in Figs. 8 and 9, in general, is a tug-of-war between 
the dominance of the symmetry producing term H r [reflecting momentum transferred from the initial conditions 
L 7 2, (P 7 2)e] an d the asymmetry producing term J?*(P 7l ) r ; more on this below. Note, like the PCS, the way in 
which the GM field acts directly on the particles in these PPP (77 — > e~e + ) processes, wherein we do not have to 
apply Eq. ip?]) manually as in these analyses, is through the presence of u>, the frame dragging velocity [cf. Eqs. ( |2^ ) 
through (|30|)], related to GM potential (see Section |II Ej ), contained in the Kerr metric [Eq. (|5j)]: The Kerr metric 
defines the paths of spacetime trajectories of particles interacting and "flowing" along geodesies, with four-momenta 
P M = (P r , Pq, L, —E), in the gravitational field (or curved spacetime) of a rotating black hole, as measured by an 
observer at infinity. 

Looking further at the different cases in Figs. 8 and 9, for the case of Figs. 8(a) and 9(a) it appears that (i 7 GM )e 
[Eq. (|47])] is somewhat balanced, as described above, in such a way that the final result is a nearly symmetric 
distribution — like that expected from a Newtonian point of view — with the ratio ~ 1.12 (defining the number of 
e~e + pairs scattered below to those scattered above the equatorial plane). In the cases of Figs. 8(b) to 8(d) the ratio 
e T ~ 1.74, ~ 2.51, ~ 2.61, respectively. Thus, in light of the general conclusion above, the distributions of Figs. 8(b) 
to 8(d) are consistent with both terms in (F' GM )e [Eq- (P^)] being significant: the first term gives the high energy PPP 
electrons a boost into opposite polar directions [cf. Fig. 1(a)]; the second term acting in dominance on the incident 
infalling photons after the PPP electrons satisfy the condition \H r L T \ > |i/*(P T ) r |, provides the asymmetry, favoring 
the positive e© direction. 

Overall, in these PPP (77 — > e~e + ) processes, we find that the intrinsic "e~e + jet" (positive e© direction) to 
counter jet (negative e© direction) asymmetry (defined by e T ) increases, first quickly, as above (cf. Fig. 8), then 
slowly, to a ratio of ~ 3 : 1, for the maximum energy that can be used for the target photon in these processes, 
E 7 2 ~ 108 GeV (i.e., before the computer simulation breaks down; see [5] for detailed construction of computer 
codes). Astrophysically, what this means is that as E 7 2 is increased, implying increasing the energy-momentum of 
the escaping PPP electrons, the luminosities in the jet (£j e t) and counter jet (Lcjet) get more one-sided, mimicking 
the effects associated with relativistic beaming, until a maximum ratio in luminosity (~ 3 : 1) is reached. We see 
this occurring when we let these PPP electrons undergo "secondary" PCS with infalling disk photons [25], allowing 
us to obtain 7-ray luminosities. It is presumed that such a ~ 3:1 ratio in jet luminosities would also exist in 
synchrotron radiation if the PPP electrons were allowed to interact with a surrounding magnetic field. Thus, it 



20 



appears that relativistic beaming near line of sight of the observer, in general, will be needed to explain observations 
if the observed Lj et /Lcjct, in a source powered by a black hole, is greater than ~ 3 for the high energy 7-ray or 
synchrotron emission. Nevertheless, importantly, these Penrose processes acting in conjunction with the GM field 
can serve the purpose of beaming the escaping particles into the polar direction, probably creating the initial jets 
of AGNs. The intrinsic magnetic field produced by the dynamo-like action associated with the swirling (vortical 
orbiting [26]) escaping Penrose produced highly relativistic plasma, about the polar axis, may be important also, at 
least in the observed synchrotron emission of the jets, and perhaps in maintaining collimation: this statement must 
be investigated further. 



In addition, when analyzing the GM radial component (F GM ) r of Eq. (47), using Eqs. (f49|), (|40|), and Fig. 2(a), 
to see what effects this component has on the PPP electrons, we get the following. For e~e + pairs scattered with 
(P T h < 0, (F GM ) r > [cf. Fig. 2(a)]. But for (P T ) e > 0, (P GM ) r > for \H*(P T ) e \ < |P e (P T )*|-as would be for 
the higher energy particles; otherwise (F GM ) r < — as would be for the lower energy particles or particles with large 
(Pf )e > 0. Since in this case, only the PPP electrons with (P T ) r > escape (on the average one of the pairs escapes, 
while the other falls into the black hole [5]), it is not clear what role the radial component of the GM force plays. 
However, looking back at the radial spectra of the PCS photons, this behavior appears to be a contributing factor 
to the progressive asymmetry of (Pp h )r seen in the escaping photons with large (Pph)e > after being scattered by 
equatorially confined target electrons [cf. Figs. 3(a)— 3(e) and 4(a)— 4(e)]; and the reestablishing the symmetry seen 
in (Pp h )r for the high energy particles [implying large (Pph)*] after being scattered by the nonequatorially confined 
target electrons [cf. Figs. 3(f)— 3(i) and 4(f)— 4(i)]. 

So, overall, for the PCS and the PPP (77 — ► e~e + ) considered here, in general, the asymmetries seen in the 
scattered particle distributions, above and below the equatorial plane, that one would, in nonrotating space (or 
pseudo flat space), expect to be symmetrical, can be attributed to the GM field acting on the incident and scattered 
particles. Specifically, it is a tug-of-war between conditions \H r P&\<\H*P r \ defined by Eq. @ that causes the 
asymmetry, which varies through a ratio of ~ 3 : 1 in the e~e + jet to counter jet for PPP (77 — ► e~e + ), and a ratio 
of ~ 5 : 1 in the photon jet to counter jet for PCS. Thus, the GM force, more or less, breaks the expected reflection 
symmetry, above and below the equatorial plane. These Penrose processes, which occur near the event horizon, allow 
us to see this general relativistic effect in action. 



IV. CONCLUSIONS 



It has been shown in this investigation that Penrose scattering processes in the ergosphere of a KBH can produce 
x-ray, 7-ray, and e~e + pair energies comparable to those observed in quasars and other AGNs. This model can work 
for any size accretion disk (thin or "thick") and any size black hole. It has been shown also that the asymmetry in 
the distribution of Penrose scattered particles, above and below the equatorial plane, can be attributed to the GM 
field — the gravitational force field produced by the angular momentum of the KBH, which causes local inertial frames 
to be dragged into the direction that the black hole is rotating (sometimes called the Lense-Thirring effect). It has 
become apparent through these Penrose processes that the GM force field lines, inside the ergosphere, close to the 
event horizon, are frame dragged into the positive azimuthal direction in the local frame of the observer at infinity, 
resulting in particle distributions appearing to break the expected reflection symmetry of the Kerr metric (which 
describes the geometry outside of a rotating black hole) above and below the equatorial plane. 

Overall, these model calculations strongly suggest that the momenta of the Penrose scattered particles naturally 
aid in the production of one-sided and two-sided polar jets, consistent with the jets associated with AGNs. Moreover, 
the variations between symmetry and asymmetry found in the scattered particle distributions, and thus in densities, 
due to the GM field acting on the space momenta of the incoming and outgoing particles, could contribute to the 
overall individual "blob-like" structures we observe in the jets of AGNs; this should, however, be investigated further. 

Importantly, we have the following conclusion from this investigation: The astrophysical conceptual fact is that 
rapidly rotating black holes give rise to Penrose processes and strong GM forces, which in turn produce high energy 
particle jets with symmetrical and asymmetrical distributions: differing by a factor as high as ~ 3 — 5, above and 
below the equatorial plane. This must be taking into account in AGN models when applying relativistic beaming 
effects [27,28,29,30]. 

In addition, this investigation suggests an answer to the long asked question concerning whether or not the asym- 
metry in the jet and counter jet distributions is intrinsic to the energy source [31]. The suggested answer, in general, 
is yes, although Doppler boosting due to relativistic beaming [32,33] when the emission axis of the jet is near the 
line of sight of the observer, may still in some cases serve as an important enhancement mechanism, to increase the 
apparent observed luminosities and energies, particularly when superluminal motion in the jet flow is apparent and/or 
the jet and counter jet luminosities differ by more than a factor ^3 — 5. Details of how these Penrose scattering 
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processes relate to AGNs, their observed luminosities, and relativistic beaming models, are discussed in other papers 
by the author [25,34]. 

Further, the model presented here to extract energy-momentum from a rotating black hole by Penrose escaping 
particles, along with the intrinsic magnetic field associated with this relativistic vortical orbiting escaping "plasma," 
which is expected to give rise to synchrotron radiation — being consistent with observations, suggests a complete, 
self-consistent model, without the necessity of the external magnetic field of the accretion disk. The attempted use of 
the external magnetic field of the accretion disk in a dominant role to extract rotational energy from a black hole, by 
having the magnetic field lines torque the horizon and/or the accretion material near the horizon (i.e., the Blandford- 
Znajek (BZ) type models; see [35] and references therein), introduces historical problems yet to be solved, namely, 
how to convert from electromagnetic energy to relativistic particle energy, such that particles are ejected along field 
lines into the polar direction; and related problems, such as how to keep the field lines from becoming tangled (or 
to be sufficiently tangled); how to generate the large magnetic field strengths needed (1) if the relativistic particles 
are to be created along the field lines and (2) to be consistent with the most luminous AGNs. Equally the causality 
problem as related to the "vacuum infinity" horizon at which upon nearing, the magnetic flux is redshifted away poses 
a problem for models assuming the magnetic field lines are "anchored" to the horizon or surrounding ergosphcric 
region; see [36]. However, the electromagnetic field of the accretion disk could, perhaps, play an important role on a 
large scale in assisting the jets of the Penrose escaping particles, in further collimation and acceleration, out to the 
observed kiloparsec distances of AGN jets. The large scale magnetic field of the accretion disk may also be important 
in transporting angular momentum outward enabling material to be accreted into the black hole [37]. Now, by the 
same token, since the GM field provides a way through these Penrose processes to transport angular momentum from 
particles inside the ergosphere, outward to infinity, it too may also provide a way to transport angular momentum 
outward, in the "static" [7] accretion disk regime (r > 2 AI, i.e., where the severe frame dragging is less): This 
statement seems worthy of an investigation. 

In addition, the Penrose process distributions described here, of escaping particles, are consistent with the disk 
returning photons [38], gravitationally focused after being emitted from radii < r ms (~ 1.2M), the marginally stable 
orbit. As measured by an observer at infinity, the PCS photons are reemitted with the highest energy photons 
scattered along geodesies concentrated in the equatorial plane. Notice, cases similar to those of Figs. 3(b), 6(b) and 
3(a), 6(a), 7(a), 7(b) can be used to explain the observations of the classical stellar/galactic black-hole candidate 
Cygnus X-l [39] and the broad Fe Ka emission line at ~ 6 keV of the bright Seyfert 1 galaxy (AGN) MCG — 

6- 30-15 [35], respectively. The case similar to that of Figs. 3(b) and 6(b) which is consistent with Cyg X-l has 
model parameters for radial infalling photons (£" p h = 3.5 keV) from a thin disk [40], that either undergo PCS by 
equatorially confined orbiting target electrons [E e ~ 0.539 MeV) or PPP (77 — > e~e + ) at r m b or r p h, respectively: 
the blueshifted energies (due to frame dragging) attained by the ~ 82% up to 90% escaping particles, returning to the 
disk to be reprocessed and/or escaping to infinity, are the following: for the PCS photons, E' h ~ 12 — 250 keV (i.e., 
up to ~ 70i?ph), with relative incoming and outgoing photon luminosities (L 7 ) ou t ~ 0.4— 130 (L 7 ); n , respectively; 
and for the relativistic PPP electrons (with E l2 ~ 5 MeV), E^p ~ 4 MeV [consistent with synchrotron radiation 
in the radio regime for B ~ 10 2 gauss, and inverse Compton scattering (of disk photons) into the hard x-rays/soft 

7- ray regime — with relative incoming and outgoing photon luminosities (£ 7 ) ou t ~ 8 — 2000 (L 7 )- m , for M ~ 30Af Q , 
at ~ 100 keV— 3 MeV, respectively] — [cf. also Figs. 8(a) and 9(a)]; see [5,25,34,41] for details. The case similar to 
those of Figs. 3(a), 6(a), 7(a), 7(b), which is consistent with MCG — 6-30-15, has model parameters for radial infalling 
photons (Epk = 1.5 keV) from a thin disk, that either undergo PCS by equatorially or nonequatorially confined 
orbiting target electrons (E e ~ 0.539 MeV or E e ~ 0.615 MeV, respectively) at r m b or PPP (77 — > e~e+) at r p h: 
the blueshifted energies (due to frame dragging) attained by the ~ 72% up to 95% escaping particles, returning to 
the disk to be reprocessed and/or escaping to infinity, are the following: for PCS photons, E' ph ~ 4.6 — 138 keV or 
E' ph ~ 6 — 169 keV for equatorially or nonequatorially confined orbiting target electrons, respectively, with relative 
photon luminosities (£ 7 ) ou t ~ 0.03 — 14(L 7 ); n or (£ T ) ut ~ 0.04 — 27(L 7 ); n , respectively; and for the relativistic 
PPP electrons (with E l2 — 6.8 MeV), E^p ~ 3.5 MeV [consistent with synchrotron radiation in the radio regime for 
B ~ 10 2 gauss, and inverse Compton scattering (of disk photons) into the hard x-ray/soft 7-ray regime — with relative 
incoming and outgoing photon luminosities (L 7 ) ou t ~ 0.07 — 20 (L 7 ); n , for M ~ 10 8 M©, at ~ 103 keV — 2 MeV, 
respectively], suggesting relatively weak, less powerful and less prominent radio jets, i.e., a radio quiet AGN, like 
a Seyfert galaxy; see [5,25,34,41] for details. Note, E l2 is assumed based on prior PCS photons with {P ph ) r < 
that satisfy conditions for the existence of a turning point at the photon orbit (discussed in Section II C ), and on 
satisfaction of the energy threshold for the particles to react (see [5]). Thus, the consistency of the results of these 
Penrose process analyses with observations strongly suggests that the appropriate way to exact energy from a KBH is 
by the "Penrose- Williams" mechanism, described here and in [5,25,34,41]: using gravity as the controlling force — to 
distinguish it from BZ-type models, which use electromagnetism as the controlling force. 

Lastly, I want to conclude with the comment, that, the physics described in this paper and the associated paper 
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of [5], can be applied to different types of high energy-momentum exchange scattering processes, i.e., processes that 
would be expected to occur in the "hot" ergosphere of a rotating black hole, the intention being to extract rotational 
energy-momentum. An immediate future investigation would be to perform the proton-proton nuclear scattering of 
neutral pions ir°, and their subsequent ir° — ► 77 decay (see [5,8,20]), inside the ergosphere of a KBH. These 7-rays are 
presumed to be important in the population of the photon orbit in the highest energy PPP (77 — ► e~e + ) processes. 
Note, similar nuclear scattering and 7r°-decay processes, however, for particles confined to the equatorial plane are 
investigated in [3]. 
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FIG. 1. Orbital parameters of massless and material particles. Magnitude of polar coordinate angular momentum Pq as 
the bound orbit crosses the equatorial plane (= Q 1 ^ 2 ) vs. the conserved nonequatorially confined orbital energy: (a) (P p h)e 
vs. E ph at the photon orbit, r ph ~ 1.074M. (b) (P e )e vs. E e , of an orbiting electron, at r mb — 1.089M (dashed curve) and 
at r ms ~ 1.2M (solid curve). Notice that, as these momenta go to zero, each of the orbital energies goes to its equatorial 
confinement value — i.e., as Q X J 2 — > 0, E pil — > 0; as Q X J' 1 — > 0, E e — > 0.539 MeV at r mb , and E e — > 0.349 MeV at r ms . 
(Note that, when the more exact value of r m b = 1.091M is used, E e — > 0.512 MeV ~ fi e -) 
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FIG. 2. The GM field force components, (a) The radial component: (/ gm )t-/p* vs. O. (b) The polar coordinate component: 
(/gm)s/p* vs - ©• The solid curves indicate the force at r m b ~ 1.089M; the dashed curves indicate the force at r p h ~ 1.074M. 
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FIG. 3. Compton scattering: scatter plots showing radial momentum components: (Pph)r vs. E' ph , of scattered escaping 
photons after 2000 events (each point represents a scattering event), at r m b — 1.089M. The various cases are defined by the 

1/2 

following parameters: E p h, initial photon energy; E e , target electron orbital energy; Q c , corresponding polar coordinate 
momentum (P e )e of the target electron; iV cs , number of photons escaping, (a) E ph = 0.511 kcV, E e = 0.539 MeV, Ql /2 = 0, 
N CB = 1707. (b) E ph = 3.5 keV, E e = 0.539 MeV, Q\' 2 = 0, N CB = 1637. (c) E ph = 0.03 MeV, E e = 0.539 MeV, 
Qo /2 = 0, N CB = 1521. (d) E ph = 0.15 MeV, E e = 0.539 MeV, Q\' 2 = 0, N CB = 1442. (e) E ph = 1 MeV, E e = 0.539 MeV, 
Qo /2 = 0, N CB = 1390. (f) E ph = 0.15 MeV, E e = 1.297 MeV, Q\' 2 = ±2.479 Mm e , N CB = 1628. (g) E ph = 0.15 MeV, 
E e = 5.927 MeV, Ql /2 = ±12.43 Mm e , N CB = 1843. (h) E ph = 0.03 MeV, E e = 5.927 MeV, Ql /2 = ±12.43 Mm e , N OB = 1935. 
(i) E ph = 0.03 MeV, E e = 11.79 MeV, Ql /2 = ±24.79 Mm e , N CB = 1971. (Note, due to a minor oversight leading to improper 
treatment in the computer simulation of the arccosine term in eq. (3.39) of [5], correct Fig. 3(i) presented here replaces Fig. 4(a) 
of [5].) 
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FIG. 4. Compton scattering: scatter plots showing the corresponding polar coordinate space momentum components: (Pp h )e 
[= (Qph) 1 ^ 2 ] vs. E' ph , of the escaping PCS photons, for the cases (a) — (i), respectively, described in Fig. 3. The units of (Pp h )e 
are Mm e (G = c = 1). (Note, due to a minor oversight leading to improper treatment in the computer simulation of the 
arccosine term in eq. (3.39) of [5], correct Figs. 4(b) and 4(c) presented here replace Figs. 7(a) and 3(c), respectively, of [5].) 
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FIG. 5. Compton scattering: scatter plots showing polar coordinate space momentum components: (-Pph)e [= (Q'ph) 1 ^ 2 ] 
vs. E' ph , of the escaping PCS photons after 2000 events (each point represents a scattering event), at r m b — 1.089M. The 

1/2 

various cases are defined by the following parameters: E p h, initial photon energy; E e , target electron orbital energy; Q e , 
corresponding polar coordinate momentum (P e )e of the target electron; N cs , number of photons escaping, (a) E p h = 0.511 keV, 
E e = 0.570 MeV, Ql' 2 = ±0.393 Mm e , N CB = 1719. (b) E pb = 0.511 keV, E e = 0.714 MeV, Ql /2 = ±0.987 Mm e , N CB = 1779. 
(c) E ph = 0.511 keV, E e = 1.297 MeV, Ql' 2 = ±2.479 Mm e , N CB = 1896. (d) E ph = 0.511 keV, E e = 5.927 MeV, 
Ql /2 = ±12.43 Mm e , iV es = 1996. 
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FIG. 6. Compton scattering: scatter plots displaying polar angles, above and below the equatorial plane: Q' ph vs. E' ph , of 
the escaping PCS photons, for the cases (a) — (i), respectively, described in Figs. 3 and 4. Note that, Op h > it/2 is below the 
equatorial plane. 
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FIG. 7. Compton scattering: scatter plots displaying polar angles, above and below the equatorial plane: Op h vs. E' ph , of the 
escaping PCS photons, for the cases (a) — (d), respectively, described in Fig. 5. Note that, Q' ph > it/2 is below the equatorial 
plane. 
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FIG. 8. Penrose pair production (77 — > e~e + ): scatter plots showing polar coordinate space momentum components: 
Pe (= Q 1 ^ 2 ) vs - E^l, of scattered escaping e _ e + pairs after 2000 events (each point represents a scattering event), at r p h. 
The various cases are defined by the following parameters: -E 7 i, the infalling photon energy; E 1 2, the target photon orbital 
energy; Q^ 2 , corresponding polar coordinate momentum (P 7 2)e of the target photon; N cs , number of e"e + pairs escaping, 
(a) E-yi = 0.03 MeV, E l2 = 13.54 MeV, Q\ /2 = ±0.393 Mm e , N as = 1850. (b) £ 7 i = 0.03 MeV, £ 7 2 = 206.7 MeV, 
7 2 2 = ±6.0 Mm e , iV cs = 1984. (c) E 7l = 0.03 MeV, £ 72 = 2.146 GeV, Q\{ 2 = ±62.28 Mm e , iV cs = 1997. (d) £ 7 i = 0.03 MeV, 
£ 7 2 = 3.893 GeV, Qlj 2 = ±113.0 Mm e , iV os = 1997. 
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FIG. 9. Penrose pair production (77 — > e e + ): scatter plots displaying polar angles, above and below the equatorial plane: 
6 T vs. E' T , of the escaping e~e + pairs, for the cases (a) — (d), respectively, described in Fig. 8. 
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